
f 




SOLUTE TRANSPORT DURING THE CYCLIC OXIDATION OF Ni-Cr-Al 
ALLOYS 


I 

0 

t 


Jj 

L3 

=3: 

Du 

r/D 

a 


CO 

C3 


CD 

ct: 


, i 


James A. Nesbitt 

Michigan Technological University 
Houghton, Michigan 


May 1982 



Prepared for 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
Lewis Research Center 
Under Grant NSG-3215 





TABLE OF CONTENTS 


Page 

I. INTRODUCTION ....... . 1 

Thermal Fatigue 2 

Oxidation Protection 2 

Philosophy and Approach 19 

Purpose 20 

II. EXPERIMENTAL RESULTS 22, 

Procedure 22 

Experimental Results 36 

Summary 67 

III. NUMERICAL MODELING RESULTS AND DISCUSSION 69 

Numerical Modeling 69 

Results and Discussion 73 

Summary 88 

IV. CONCLUSIONS 90 

APPENDICES 91 

Appendix A: Technique of Quantitative Electron Microscopy 91 

Appendix B: Effective NiCrAl Phase Diagram at 1200°C ... 95 

Appendix C: Diffusion Coefficients . 100 

Appendix D: The Numerical Model 118 

REFERENCES 134 


i 



1 


Solute Transport During the Cyclic Oxidation of Ni -Cr-Al A1 1 oys 

by James A. Nesbitt 

Michigan Technological University 
Houghton, Michigan 

INTRODUCTION 

Most superalloys require protective surface coatings to achieve practical 
component lives in high-temperature, oxidizing, and corrosive environments. 
Continued advances in the high- temperature mechanical properties of Ni and Co 
based superalloys have increased their possible operating temperature to 1200°C 
(1). Ni or Co alone do not form oxides capable of surface protection at the 
elevated temperatures (2). In the past, many superalloys contained sufficient 
quantities of Cr to provide protective Cr202 oxide scales whereas the recent 
advances in the mechanical properties have come mainly at the expense of the 
Cr content (3). The result has been an increase in the structural stability 
of superalloys while the surface stability has decreased (2). 

Protective coatings are generally of two types, aluminide diffusion 
coatings and overlay coatings. Currently, the m.ost promising overlay coatings 
are alloys in the NiCrAl and CoCrAl systems (1). Both coating types provide 
high- temperature oxidation protection by the selective oxidation of A1 . 

3+ 

“”^"*2^3’ high-temperature oxide phase, allows minimal transport of A1 

2 _ 

and 0 ions, thereby providing sufficient oxide protection between the metal 
and the environment. In general, the most important requirements for protec- 
tive coatings for gas turbine superalloys are (1): 

1. resistance to oxidation accompanied by thermal cycling, 

2. resistance to thermal fatigue cracking, and 

3. resistance to hot corrosion (for specific applications and environ- 
ments). 

The bulk of this report will deal with resistance to oxi dation accompanied 
by thermal cycl ing. The resistance to thermal fatigue cracking will be briefly 
discussed below. The resistance to hot corrosion will not be considered. 
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Thermal Fatigue 

Thermal fatigue results from mechanical strains induced in superalloy 
components undergoing thermal cycling. Many superalloy appl icationSj such as 
jet engine turbine components, entail thermal cycling during normal use. The 
danger of thermal fatigue is component cracking, especially in brittle surface 
coatings (4). The strains are typically great enough to propagate cracks into 
the superalloy component (1) resulting in catastrophic failure. Hence, some 
ductility is required in protective coatings to resist thermal fatigue cracking. 

Two-phase overlay coatings are capable of greater ductility than aluminide 
diffusion coatings. The surface phase of aluminide coatings is based on the 
intermetal 1 ic NiAl (CoAl) phase. It is unfortunate that the NiAl phase is 
quite brittle below approximately 700 - 900°C depending on the exact A1 concen- 
tration (3,4). For many advanced airfoil designs, thermal strains can peak at 
low temperatures, and brittle coatings can prove to be inadequate for some 
applications (1). In comparison, two-phase NiCrAl (CoCrAl) overlay coatings 
combine the oxidation resistance of the brittle NiAl (CoAl) phase with the 
ductility of the Ni (Co) solid solution phase. Naturally, the ductility of 
the two-phase structure is dependent on the amount of the NiAl (CoAl) phase 
present. Therefore, two-phase overlay coatings sacrifice A1 content (less 
NiAl phase) for the required greater ductility (1). 

Oxidation Protection 

Oxide protection by overlay coatings depends on the ability of the coating 
to supply continuously the A1 necessary for growth. When, 

the scale is adherent and intact, the slow ionic transport in the 

scale is the rate-controlling step (5), thereby fixing the overall rate of A1 
consumption. The ionic transport, or diffusion, in the oxide scale is depen- 
dent on the diffusivities of the ions, the interface concentrations at the 
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gas/oxide and oxide/metal boundaries, and the thickness of the oxide scale (6). 

As the scale thickens, the ions must travel greater distances, and the rate of 
A1 consumption decreases. 

Control of the oxidation rate by diffusion in the scale allows definition 
of a specific flux or delivery rate of A1 , viz., the requisite flux. The 
requisite flux is the specific delivery rate of A1 from the metal to the oxide 
required for exclusive formation and continued growth of an Al^O^ oxide scale. ^ 

As so defined, the requisite flux is equal to the rate of A1 consumption (7) 
as controlled by diffusion in the oxide scale and is therefore dependent on 
the oxide scale thickness. When the oxide scale is thin, the requisite flux 
will be great; as the oxide scale thickens, the requisite flux will decrease. 
Oxide protection by the overlay coating is afforded when the flux of A1 from 
the coating substrate to the oxide scale is equal to the requisite flux. 

Overlay coatings when bonded to superalloy substrates undergo loss of A1 
by interdiffusion with the substrate (1,9). A1 diffusion into the superalloy 
substrate effectively reduces the A1 content of the coating. The degree of 
interdiffusion will depend on the diffusivity of A1 in the superalloy substrate 
(9,10). If the A1 diffusivity is high, more A1 can be lost to the substrate 
than is consumed in the formation and growth of the Al^O^ oxide scale (11). 

As the A1 reserves in the coating are depleted, the coating will be unable to 
supply the requisite flux to the oxide/coating interface, and the coating may 
cease to provide adequate protection (10,11). 

The diffusional processes involved with the oxidation of coated superalloys 
may be divided into two related studies: the diffusional transport of A1 to the 
oxide/coating interface, and the interdiffusion of A1 between the coating and 

The requisite flux as defined above is equivalent to the ionic flux 
designated S as used by Wagner (8) and the "flux of A1 consumed due to oxidation, 

discussed by Pettit (7). 

oxide 
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superalloy substrate. This thesis will deal mainly with the former study, viz., 
the diffusion of A1 to the oxide/metal interface, specifically in two-phase 
NiCrAl alloys. 

Thermal cycling of NiCrAl alloys results in cracking and spalling of the 
protective Al^O^ oxide scale (12). The oxide cracking and spalling is a con- 
sequence of mismatch between the coefficients of thermal expansion of the oxide 
and the metal substrate (13,14). Since the requisite flux is greater for thin 
oxide scales, the requisite flux associated with oxide spallation will be 
greater than if the oxide scale had remained adherent and intact. As a result, 
the rate of A1 consumption will also be greater when cracking and spalling of 
the oxide scale occurs (3,15). The danger of a higher rate of A1 consumption 
is a more rapid depletion of the A1 reservoir in the two-phase region. When 
the alloy is unable to supply the requisite A1 flux to the oxide/metal inter- 
face, the poorly protective oxides of Ni, Co, or Cr will form (7,15). This 
failure of the protective oxide scale is commonly referred to as breakaway 
oxidation (6) and occurs in many alloys where a reactive metal (as Cr, Si, or 
Al) is selectively oxidized to form a protective oxide scale (16,17). 

The ability of the alloy to supply the requisite flux is dependent on the 
solute diffusivity and the concentration gradients near the oxide/metal inter- 
face. The effect of the diffusivity and the bulk alloy composition (with 
regard to the phase diagram) on the ability to supply the requisite flux in 
two-phase alloys can be illustrated most simply by reference to three theoret- 
ical binary alloys. Assume that all three alloys form the same solute oxide 
at the same rate so that the amount of solute consumed at any time will be the 
same. For simplicity, imagine the oxide to remain at a constant thickness 

^In the remainder of this report, metal, metal substrat¥, and alloy will 
refer to the coating alloy unless otherwise specified. 
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(as results from oxide volatization or spallation occurring at the same rate 
as oxide growth). Therefore, the requisite flux for each alloy is constant 
and equal. 

For the first alloy (alloy 1), consider a slab of a y + B> A-B alloy where 
g is the B rich solute phase, and y is the A rich solvent phase. If a slab of 
alloy 1 is oxidized, B will diffuse to the oxide/metal interface where it is 
oxidized, and a y layer depleted of the g phase will grow inward from the metal 
surface. A schematic drawing of the B concentration/distance profile for two 
oxidation times, t < t* and t=t*, is shown in Figures la and lb. The initial 
solute concentration, C^, and the maximum B solubility in the y phase, C^g, 
are shown in a possible A-B phase diagram in Figure le. For simplicity only, 
the B concentration profiles in the near-surface y layer have been drawn as 
straight lines. The total amount of B consumed after t* hours can be deter- 
mined as the difference between the initial B concentration profile and the B 
concentration profile at t = t*, indicated by the cross-hatched area in Figure 
lb. 

The B concentration gradients at t<t* and t = t* (Figures la and lb) are 
equal because at both times the alloy is supplying the constant requisite flux. 
Up to the time t=t*, as the width of the layer increased, the B concentration 
at the oxide/metal interface could decrease, maintaining a constant B concen- 
tration gradient and supplying the constant requisite flux. But at t=t*, the 
B concentration at the oxide/metal interface has decreased to zero (Figure lb). 
As the y-layer width continued to increase (t>t*), the B concentration gradient 
must decrease, alloy 1 would no longer be capable of supplying the requisite 
flux, and the B oxide scale would fail; i.e., oxides of the solvent A would 
begin to form.^ 

^Basic thermodynamic considerations show that the A oxides begin to form 
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Figure 1. Oxidation of y+S Binary Alloys 
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Increasing the diffusivity enables alloy 1 to continue to supply the 
requisite flux and provide a protective B oxide scale for a longer time (t>t*) 
Doubl i hg the diffusivity of alloy 1 would allow the B concentration gradient 
to be reduced by a half and yet still supply the requisite flux. At t = t*, 
the y- layer width of alloy 1 with twice the diffusivity would be slightly 
greater than that shown in Figure lb, but the B concentration at the oxide/ 
metal interface would be significantly greater than zero (due to the lower 
concentration gradient). Hence, increasing the diffusivity increases the 
ability of an alloy to supply the requisite flux. The drawback to an increased 
diffusivity is the rapid loss of solute from a coating to the metal substrate, 
more rapidly depleting the coating and possibly resulting in an earlier oxide 
scale fai 1 ure (t < t*) . 

The effect of the initial solute concentration on the solute flux can be 
illustrated by considering a second y+8 alloy (alloy 2) with a greater B 
content (Cq > Cq; D, same as alloy 1) oxidized for t* hours. The resulting 
B concentration/distance profile is shown in Figure Ic. Since the amount of B 
consumed is controlled by the requisite flux of the oxide scale, both alloys 1 
and 2 have equal concentration gradients in the y layer and also have lost the 
same amount of B at t = t* (the cross-hatched areas in Figures lb and Ic are 
equal). Since ^ , the y- layer width in alloy 2, is less than the y-1 ayer 
width in alloy 1, the B concentration at the oxide/metal interface in alloy 2 
is greater than zero. Hence at t " t* , alloy 2 does not exhibit breakaway 
oxidation but is still able to supply the requisite flux to the oxi de . There- 
fore, increasi ng the sol ute concentration increases the protecti ve 1 ife of 

before the B concentration at the oxi de/metal interface is equal to zero. But, 
for a solute with a very high affinity for oxygen (as in the case of A1 in Ni 
or Co) , the same thermodynamic considerations show that the solute concentra- 
tion at the interface is many orders of magnitude less than the solvent concen- 
trati on at the same interface (i .e. , 1 ppm A1 in Ni ) (7) and may therefore be 
approximated as zero. 
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oxide scales on binary alloys. 

A third alloy (alloy 3) with the same B concentration but a greater B 
solubility in the y phase (C*^ > Cq, D same as alloy 1) can be used to 
illustrate the beneficial effect of the greater B solubility. After oxidation 
for t* hours, alloy 3 is still able to supply the requisite flux. The B con- 
centration gradients in the y layer of both alloys are equal, but the B con- 
centration at the oxide/metal interface in alloy 3 is greater than zero 
(Figure Id), The detrimental effect of the opposite change from that illus- 
trated above for Cq, or D could be similarly reasoned. In a ternary alloy 
system, the initial solute concentration (Cq), the maximum solute concentration 
in the solvent phase (C^g), and the diffusivity (D) may all vary at the same 
temperature. The detrimental effect of two of the above parameters could over- 
come the beneficial effect of the third. For example, it is possible to 
envision a ternary alloy with a greater initial solute concentration failing 
before another ternary alloy with a lower initial solute concentration. 

Oxidation studies can be arbitrarily divided into tv;o overlapping cate- 
gories, isothermal oxidation and oxidation accompanied by thermal cycling. 

The main emphasis of this report concerns the latter category dealing with the 
A1 diffusion in two-phase NiCrAl alloys during oxidation accompanied by thermal 
cycling. Since isothermal oxidation studies provide insight to the oxidation 
behavior observed during the more complex cyclic oxidation,^ the pertinent 
isothermal oxidation studies will be briefly reviewed before considering cyclic 
oxidation. 


^Cyclic oxidation will hereafter refer to oxidation accompanied by 
thermal cycl ing. 
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Isothermal Oxidation 

Experimental . High temperature isothermal oxidation of (Ni ,COjFe)CrAl 
alloys results in the selective oxidation of A1 (5,18,19). Alloys containing 
greater than 4-6 wt.%Al form external AI2O2 oxide scales which generally 
thicken parabolically with time (5,18). Since the oxide growth is controlled 
by diffusion in the oxide itself, there is little difference between the growth 
rates of external AI2O2 oxide scales on various base alloys (Ni,Co, or Fe) (18) 
or alloy composition. The differences appear negligible in comparison with 
the growth rates of other possible oxides (NiO, CoO, or Cr202) (3,5,18,20). 

After short oxidation periods, transient oxides of Ni, Co, Fe, or Cr 
typically appear on the alloy surface (5,21-23). The greater the A1 concen- 
tration of the alloy, the quicker the replacement of transient oxides by a 
continuous A1202 oxide scale. For alloys with low (less than 4 wt.%) A1 con- 
centrations, the less-protective oxides may never be fully replaced by a con- 
tinuous AI2O2 oxide scale. When the less-protective oxides remain on the 
alloy surface, oxygen diffusing into the alloy internally oxidizes the A1 with 
the resultant growth of a subscale of A1202’ precipitates below the sample 
surface (5,18). 

The transition from transient oxide formation to a continuous AI2O3 oxide 
scale and the absence of such a transition is directly related to the diffusion 
in the near-surface region. Alloys with high A1 concentrations are rapidly 
able to develop large A1 concentration gradients and supply the relatively 
high A1 flux necessary to establish an external scale of AI2O2. In comparison, 
alloys with less than 4-6 wt.%Al cannot develop the flux necessary for forma- 
tion of an external layer of AI2O2. Internal oxidation of A1 occurs below 
the less-protective oxides (5,18). The denser the AI2O3 precipitation in 
the subscale layer, the more tortuous the path the incoming oxygen must travel. 
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and the greater the chance of formation of a continuous layer below 

the subscale region. If a significantly large volume fraction of Al^O^ pre- 
cipitates is deposited at the subscale/metal interface, the Al^O^ precipitates 
may "link up," forming a continuous AI 2 O 2 layer (5,24). If a continuous 
layer of AI 2 O 2 is not formed, continued internal oxidation essentially depletes 
the alloy of free A1 . The resulting oxidation rate is relatively high, being 
associated with the less protective oxides Cr^O^, NiO, CoO (5,18), etc. 

Oxide maps, experimentally determined, indicate regions of specific oxide 
control on isotherms of the ternary phase diagram. Oxide maps for the NiCrAl 
system at 1200°C (5) showing regions of Al^O^, Cr^O^, NiO, and spinel oxide 
formation areshown in Figure 2. 

Small additions of reactive metals (as Y and Zr) greatly increase the 

A1 0„ oxide scale adherence of (Ni ,Co,Fe)CrAl alloys (19,25-27). Most proposed 
2 d 

coating alloys for high-temperature use, therefore, contain various amounts 
(< 1 wt.%) of these metal additions. Various proposed mechanisms to account 
for the increased adherence have been discussed elsewhere (19,25,28-31). 
Numerous isothermal oxidation tests have shown that the reactive metal addi- 
tions can significantly affect the oxidation rate (19,25,27). Oxide stringers 
growing from the oxide scale into the metal substrate are observed on oxidized 
alloys containing reactive metal additions. The origin or cause of the oxide 
stringers is unclear for Zr-containing alloys where no enrichment of Zr is 
observed (27). The oxide stringers observed in Y-containing alloys are attrib- 
uted to the preferential oxidation of grain boundary yttrides (19,25,28,29). 

Theoretical . A theoretical solute concentration/distance profile for 
oxidized alloys can be predicted by an appropriate solution to Fick's second 
law subject to two boundary conditions. One boundary condition may be the 
initial solute concentration in single-phase alloys or an interface 
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50 at .®/o 
Al 



Figure 2 Oxide maps for the isothermal and cyclic oxidation of NiCrAl 
alloys at 1200®C. Solid lines indicate limits for isothermal 
oxidation, dashed lines indicate limits for cyclic oxidation. 
References 5 and 12. 
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concentration when a second-phase is involved. The second boundary condition 
at the oxide/metal interface may be either a solute concentration or a flux 
condition. For predicting the solute concentration/distance profiles after 
oxidation of y + S alloys, sufficient boundary conditions would be the initial 
solute concentration (Cg(t)) and the maximum solute solubility (C^g(t)) at the 
y/y + 6 interface, and either the requisite flux (J(t)) or the solute concentra- 
tion (C^(t)) at the oxide/metal interface (see Figure 1). Wagner (8) has shown 
that the solute concentration at the oxide/metal interface is not sufficiently 
determined by thermodynamic considerations alone, but may be found by equating 
the solute flux in the alloy with the ionic solute flux in the oxide. Wagner 
has also shown that when the oxidation rate is controlled by diffusion in the 
oxide, the solute concentration in the alloy at the oxide/metal interface is 
non-zero. Ku (32) has shown that if the oxide scale grows parabol ically , the 
oxide/metal interface concentration is independent of time."^ The time- 
independent solute concentration at the oxide/metal interface during parabolic 
oxidation can be compared with the time-dependent solute concentration at the 
same interface during oxidation when the oxide scale thickness and requisite 
flux are constant, as shown in Figure 1. Bastow et al. (33) have predicted 
time-dependent interface concentrations during the transient oxidation period 
when non-parabolic rate laws apply. 

Various researchers have applied a diffusional analysis to predict the 
concentration/distance profiles of binary alloys exhibiting parabolic oxidation 

^Parabolic oxidation, parabolic oxidation behavior, and parabolic oxida- 
tion kinetics will refer hereafter to oxidation which results in a parabolic 
growth of the oxide scale (^^oxide ” ^ parabolic decrease in the requisite 

flux (J = ko//t) , a time-independent solute concentration at the oxide/metal 
req c. 

interface (Cs = kT), and for y+B type alloys, a parabolic growth of the near- 
surface Y layer (AX^ = k^/t) . 
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kinetics. FeCr (33-36) and NiCr (37,38) alloys have received the bulk of the 
research. A paucity of the necessary thermodynamic and diffusion data has 
precluded a diffusional analysis of the oxidation of ternary alloys. 

Cyclic Oxidation 

Experimental . Cyclic oxidation studies have concentrated on two areas: 
investigating the various parameters affecting oxide spallation, and deter- 
mining the effect of oxide spallation on the overall oxidation and metal con- 
sumption rate. Early studies (20) have shown Al202/NiAl20^ oxides or certain 
Cr 202 oxides to possess the best cyclic oxidation resistance, i.e., to be the 
most adherent scales. However, Cr20^ oxide scales are not considered practi- 
cal due to transformation of Cr 202 to CrO^* highly volatile at temperatures in 
excess of 1100°C (2,3). The excessive volatization which occurs with high gas 
velocities or low pressures could lead to early breakaway oxidation (20). 
Alumina-forming alloys have been chosen as the best candidates for cyclic 
applications at temperatures in excess of 1100°C. An oxide map indicating the 
regions of AI 2 O 2 , Cr202» NiO, and spinel oxide formation during cyclic oxida- 
tion of NiCrAl alloys at 1200°C (12) is shown in Figure 2. 

Numerous investigations have shown oxide spallation to be affected signifi 
cantly by several parameters. Those parameters include cooling rates (39), 
quench temperature (40), oxide dispersions in the alloy (41-43), and the cycle 
frequency (the duration of the heating and cooling cycle) (17,44). Small addi- 
tions of Pt to the alloy, or the presence of Pt in the furnace, have also been 
observed to increase the oxide adherence (45-47). Probably the most signifi- 
cant and beneficial effect is that produced by small additions of reactive 
metals. Zr was found to improve significantly the scale adherence of alumina- 
forming NiCrAl (12,48) and CoCrAl (44) alloys, and chromia-forming FeCr (41) 
and Cr stainless steels (49), Y has also been found to increase the scale 



15 


adherence of A1202 on NiCrAl (19,25,30), CoCrAl (19,50), and FeCrAI (26,28,29, 
51) alloys. In addition. Sc (28), Th (25), and Hf (50) have been shown to 
exhibit a beneficial effect on scale adherence. 

Evaluation of various alloys and active metal additions is typically 
accomplished by comparison of sample weight change measurements'*' made during 
cyclic oxidation. During the high- temperature portion of a thermal cycle, the 
sample weight increases parabolically with time as the oxide scale thickens. 
Upon cooling, patches of oxide spall non-uni formly, decreasing the weight of 
the sample. Depending on the scale adhesion, the overall weight change will 
eventually become negative. When the oxide growth is approximately equal to 
the oxide loss during each cycle, the weight change will decrease linearly with 
increasing oxidation cycles. Alloys undergoing cyclic oxidation accompanied by 
scale spallation are said to exhibit "paral inear-type" (6) oxidation behavior, 
descriptive of the weight change when plotted as a function of time as shown in 
Figure 3 (alloy 5A). (Paral inear- type oxidation behavior is in contrast to the 
parabolic weight increases observed during isothermal oxidation.) When break- 
away oxidation occurs, the less protective oxides which form tend to spall at 
greatly accelerated rates causing an even more rapid decrease in the weight 
change (17), also shown in Figure 3 (alloy 14). 

Theoreti cal . Two modeling techniques are typically employed to determine 
the weight of metal consumed during cyclic oxidation. The two techniques are: 

1. Macroscopic spalling models, and 

2. Microscopic diffusion models. 

Macroscopic spalling models predict the total weight of oxide formed, the 
weight of the oxide which has spalled, and the weight of oxide remaining 

. 4 . 

For comparison of various sample sizes, the weight change is typically 
divided by the total sample surface area (mg/cm^). 
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Figure 3 



Comparison of the cyclic oxidation behcivior of CoCrAl 
alloy 14 (Co-18. 93at,%Cr-23.25at.%Al) and NiCrAl alloy 
5A (Ni-14.35at.%Cr-23.65at.%Al) at 1200^ C for 2000 0.1 
hour exposure cycles. 5A data are paralinear in nature 
with scale of Al 2®3 aluminate plus ZrO^ detected. 

Alloy 14 data breaR at close to 115 hour due to CoO 
buildup and subsequent spalling (Peference 44) . 
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intact on the satrple surface. Input to most spalling models consists of an 
oxide growth constant and various spalling parameters to govern the amount of 
oxide which spalls each cycle. Various models have been proposed for AI 2 O 2 
spallation from NiAl (52) and NiCrAl (53) alloys, and for the analogous case 
of oxide volatization from Cr^O^-forming alloys (54). 

Microscopic diffusion models simulate the solute diffusion to the oxide/ 
metal interface. The instantaneous solute flux at the oxide/metal interface 
is either known, and input to the model as a boundary condition, or can be 
determined from the solute concentration/distance profiles predicted by the 
model. From the instantaneous solute flux, and the predicted concentration/ 
distance profiles, a diffusion model has the inherent capability to predict 
both the weight of solute consumed as oxide and the time to the onset of 
breakaway oxidation. 

The solute concentration/distance profiles during cyclic oxidation are 
described by a solution to Pick's second law subject to tv/o boundary condi- 
tions, the same as for the case of isothermal oxidation. The inner boundary 
condition is also the same as that for isothermal oxidation, viz., either an 
initial solute concentration or interface concentration for one-phase or multi- 
phase alloys, respectively. The significant difference in predicting the 
concentration/distance profiles during cyclic oxidation and isothermal oxida- 
tion is the nature of the boundary condition at the oxide/metal interface. 

The boundary condition at the oxide/metal interface during cyclic oxida- 
tion is critically dependent on the extent of oxide spallation. As previously 
stated, the requisite solute flux is inversely related to the oxide scale 
thickness. The thinner the oxide scale, the greater the requisite solute flux. 
During cyclic oxidation accompanied by oxide spallation, the average oxide 
thickness will always be less than if the scale had been fully intact. As a 
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result, the requisite flux*^ of alloys undergoing oxide spallation during 
thermal cycling will always be greater than the requisite flux associated with 
intact, adherent scales. The greater the extent of oxide spallation, the 
greater the requisite flux. 

An alloy must supply the requisite flux if the oxide scale is to remain 
protective during cyclic oxidation. For the case of parabolic oxidation of 
a single-phase binary alloy, Wagner (8) and others (33,35-37) have shown that 
balancing the ionic flux in the oxide with the solute flux in the alloy 
results in a solute concentration at the oxide/metal interface less than the 
initial bulk solute concentration of the alloy. Therefore, a solute concentra- 
tion gradient is established which supplies solute from the metal substrate to 
the oxide. The solute concentration at the oxide/metal interface remains 
constant as long as the oxide scale thickens parabol ical ly (32). But, if 
during isothermal oxidation all or part of the oxide scale were to spall off, 
rebalancing the ionic oxide flux and the solute alloy flux would result in a 
further decrease in the solute concentration at the oxide/metal interface. 
Whittle (55) has derived an approximate equation describing the solute concen- 
tration in a single-phase binary alloy after total spallation of the oxide 
scale. The solute concentration at the oxide/metal interface, as calculated- 
by Whittle, is always less than the interfacial solute concentration before 
the oxide scale was lost. 

Cyclic oxidation results in repeated oxide growth followed by total, or 
partial, oxide spallation. Repeated application of the previous argument 
would indicate a discontinuous decrease in the solute concentration at the 
oxide/metal interface during cyclic oxidation accompanied by scale spallation. 

Due to non-uniform oxide spallation, the requisite flux during cyclic 
oxidation will always be considered as an average flux evaluated across a 
sanple surface. 
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The extent of the decrease in the interfacial solute concentration would be 
expected to depend on whether all, or only a portion, of the oxide scale 
spalled. Similarly, the previous statements could be extended to predict a 
decreasing interfacial solute concentration during the cyclic oxidation of 
two-phase alloys. 

The solute concentration at the oxide/metal interface must not be depleted 
below that for thermodynamic stability of the solute oxide if the alloy is to 
continue to provide a protective oxide scale (55). Stated otherwise, break- 
away oxidation occurs during cyclic oxidation when the interfacial solute 
concentration has decreased to approximately zero. Time-dependent interfacial 
solute concentrations resulting in breakaway oxidation were shown schematically 
in Figures la and lb. 

Some researchers have predicted time-dependent oxide/metal interface 
concentrations by assuming various oxide growth rate laws. Most of the pre- 
dictions involve the non-parabolic oxidation rates associated with the initial 
transient oxidation of binary alloys (56), Bastow et al . (33) have extended 
the analysis to include the prediction of breakaway oxidation in a CoCr alloy 
exhibiting paralinear oxidation behavior. The decrease in the Cr concentra- 
tion at the oxide/metal interface was predicted as a result of the non- 
parabolic requisite flux associated with Cr 202 vol atization . 

Philosophy and Approach 

The cyclic oxidation of y+B , NiCrAlZr(Y) alloys results in the selec- 
tive oxidation of Al to form a protective Al^O^ oxide scale. A y-phase sub- 
surface layer thickens with time as Al is transported to the growing oxide. 
Active metal additions affect the rate of oxide spallation which governs the 
requisite Al flux for the sole formation and growth of AI 2 O 2 . The 
requisite Al flux fixes the Al transport in the y-phase layer which controls 
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the growth of the layer. The ability of the alloy to supply the requisite 
flux is dependent on the diffusivity of the particular alloy, the maximum 
possible A1 gradient in the y layer, and the effect of the Cr concentration 
profile on the diffusion of A1 . The maximum possible A1 gradient is a function 
of the alloy composition with respect to the NiCrAl ternary phase diagram. 
Hence, varying the Y or Zr content of the alloys affects the rate of oxide 
spallation, which in turn affects the requisite A1 flux, and therefore the A1 
transport and growth rate of the y-phase subsurface layer. 

Purpose 

Verification 

The first goal of this study was to verify the effect of various requisite 
flux conditions on the A1 transport in y+g, NiCrAlZr(Y) alloys undergoing 
cyclic oxidation. The y+g alloys were chosen so as to allow examination of 
the effect of various concentrations, various maximum solubilities of A1 in 
the y phase, and various diffusivities on the diffusional transport of A1 and 
on the ability of each alloy to avoid breakaway oxidation. Varying the Zr or 
Y content of the alloys was used to induce differing oxide spallation char- 
acteristics. The effect of the above parameters was examined by optical 
microscopy and quantitative electron microscopy of the alloys after various 
cyclic oxidation exposures. X-ray diffraction of the oxide scale retained on 
the sample surface was used to identify the oxide phases. 

Prediction 

The second goal of this study was development of a numerical model to 
simulate the A1 transport during the cyclic oxidation of y+g, NiCrAlZr(Y) 
alloys. Diffusion models are capable of predicting concentration/distance 
profiles and phase layer widths (in multiphase applications). Diffusion 


V» 
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models typically require the input of boundary conditions and diffusivities 
specific to each alloy. The accuracy and validity of the predictions by the 
diffusion model would be ascertained by a direct comparison with the results 
from the experimental analysis of the cyclically oxidized alloys. 
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EXPERIMENTAL RESULTS 


Procedurg 

Three alloy compositions were chosen for examination of the A1 transport 
during the cyclic oxidation of y+3 NiCrAl alloys at 1200°C. Two alloy compo- 
sitions contained approximately the same A1 concentration but different Ni and 
Cr concentrations, whereas the third alloy composition contained a greater A1 
concentration but lower Ni and Cr concentrations. Different Zr contents were 
added to each of the alloys to effect various oxide spalling characteristics. 
The Zr was added by one of two techniques. One technique consisted of induc- 
tion melting the metallic elemental constituents Ni , Cr, A1 , and Zr in 
alumina crucibles. In the second technique, metallic Ni , Cr, and A1 were 
induction melted in zirconia crucibles which permitted Zr to be leached 
from the crucible. The longer the melt was held in the zirconia crucible, the 
greater the Zr content of the resulting alloy. One alloy composition was 
alloyed with Zr by both techniques to examine the effect of the Zr content and 

-t* 

alloying technique on the oxidation behavior. The same alloy composition was 
also alloyed with Y to allow a specific comparison of the effect of Y and Zr 
on the diffusional transport of A1 during cyclic oxidation. Following the 
induction melting, all the alloy melts were poured into zirconia leaf 
molds. Very little Zr was leached from these molds (57). The compositions of 
the alloys and the method of alloying are listed in Table I. 

The alloys were fabricated and oxidized at the Lewis Research Center of 
the National Aeronautics and Space Administration. The alloys were oxidized 
in the cast state without previous heat treatments. Isothermal oxidation was 

^Different alloying techniques and different alloy additions (Zr or Y) to 
the same alloy composition resulted in variations in the A1 and Cr concentra- 
tion (of 1-2 at.%) between the alloys. 
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TABLE I 

ALLOY COMPOSITIONS 


Conposition (at.%) 


Alloying 


Alloy Designation 

Ni 

Cr 


Zr(Y) 

Technique 

Ni-13Cr-25Al-0.05Zr 

61.95 

13.26 

24.73 

0.054 

a 

Ni-14Cr-25Al-0 . 07Zr 

61.92 

14.35 

23.66 

0.066 

b 

Ni-14Cr-25Al-0. 12Y 

61.16 

13.75 

24.97 

0.12Y 

a 

Ni-19Cr-24Al-0.04Zr 

57.37 

19.00 

23.59 

0.038 

. b 

Ni-17Cr-28Al-0. 18Zr 

54.41 

17.06 

28.35 

0.18 

b 


a. Alloyed with Zr or Y by induction melting of the elemental, metallic 
constituents in an alumina crucible. 

b. Alloyed with Zr by induction melting metallic Ni, Cr, and A1 in a 
zirconia crucible allowing Zr to be leached from the crucible. 
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performed at 1200°C in still air. Cyclic oxidation testing consisted of a 
one-hour exposure at 1200°C in still air followed by a minimum cooling period 
of twenty minutes at approximately room temperature, also in air. The oxida- 

■ 4 * 

tion/cooling cycle was repeated various times for each of the alloys. The 
change in weight per unit sample area was measured during both the isothermal 
and cyclic oxidation tests. Paral inear weight change curves of the four Zr- 
containing alloys indicate a significant range in the oxide spalling character- 
istics between the alloys. The weight change curves are shown in Figure 4. 

At longer times, the Ni-13Cr-25Al-0.05Zr and Ni-17Cr-28Al -0. 18Zr alloys exhibit 
the most rapid weight loss , indicating poor scale adherence. In comparison, the 
Ni-19Cr-24Al-0.04Zr and the Ni-14Cr-24Al-0.07Zr alloys exhibit a much lower 
weight loss, indicating better scale adherence. Further details of alloy 
fabrication, the oxidation equipment, and testing procedures are given else- 
where (12,27,48). 

The three alloys containing 24 - 25 at.%Al initially possess very similar 
Y + 8 microstructures with 0.6-0.75 volume fraction of the g phase present. 

The Ni-17Cr-28Al-0. 18Zr alloy initially contains approximately 0.96 volume 
fraction of the g phase. (The volume fractions of the g phase decrease after 
long cyclic oxidation exposures. The microstructural changes during cyclic 
oxidation will be discussed shortly.) The microstructures of the four alloys 
and the alloy position in the y+g field of the NiCrAl ternary diagram at 
1200°C are shown in Figure 5. The dendritic g phase (dark phase in the micro- 
graphs of Figure 5) has the NiAl (or CsCl ) type crystal structure. The y phase 
(light phase in the micrographs of Figure 5) is a Ni solid solution with an FCC 

**"The cyclic oxidation time referred to in this report will only account 
for the high-temperature period of each cycle. Hence, cyclically oxidized 
for 800 hours refers to 800 cycles, each consisting of one hour at 1200°C 
followed by a minimum of 20 minutes at room temperature. 



AW-Weight Change (mg/cm^) 



One-Hour Oxidation Cycles 


Figvire 4. Paralinear weight change curves for the Zr-containing 

alloys during cyclic oxidation. Parabolic weight gains 
during isothermal oxidation of the Ni-14Cr-24Al-0.07Zr 
alloy is also shown. 




Ni - 19.0 Cr - 23.6 Al-.042r Ni -13.8 Cr -25.0AI - .1 2Y 


Pigrare 5. Microstructures of the foxir, y+g NiCrAl alloys 
examined in this study. 



27 


type crystal structure. No enrichment of Zr was detected by optical metallo- 
graphy, x-ray mapping on the electron microscope, or by x-ray diffraction of 
the three alloys containing Zr. The Y-containing alloy contained a low volume 
fraction of a Ni and Y enriched phase at the boundaries between the y and 3 
phases, most probably Ni yttrides (19,25,28,29). Hence, the solubility limit 
of Y was exceeded but the solubility limit of Zr was not, in accordance with 
previous studies (19,25,27,30). 

The oxide phases retained on the surface of each oxidized alloy were 
determined by x-ray diffraction. The intensities of the oxide and metal sub- 
strate peaks were ranked to indicate the proportion of each of the oxide 
phases present. The oxide and metal phases detected are listed and ranked in 
Table II. The Ni-14Cr-24Al-0.07Zr and Ni-19Cr-24Al-0.04Zr alloys show the 
strongest AI 2 O 2 presence through 800 hours. Both alloys exhibited good oxide 
spalling characteristics evident from the weight change data in Figure 4. The 
Ni-13Cr-25Al-0.05Zr and Ni-17Cr-28Al -0. 18Zr alloys show considerable spinel 
formation after 500 hours of cyclic oxidation exposure. Both alloys exhibit 
poor oxide spalling characteristics, also evident from the weight change data 
in Figure 4. Continuous formation of AI 2 O 2 results in the best oxide spalling 
characteristics. 

The oxide morphologies of the four alloys are as varied as the paral inear 
weight change data of Figure 4. The oxide morphologies differ between alloys 
and with increasing oxidation time. Examples of the various morphologies are 
shown in Figure 6. The oxide morphologies range from relatively flat (Figure 
6a) to rough and porous (Figures 6b and 6g). Oxide spallation from a sample 
surface does not occur uniformly but from random, discrete sections on the 
sample surface. The sample surface of the Ni-14Cr-25Al-0. 12Y alloy, shown in 
Figure 7, illustrates the random oxide spallation. The light areas in 



TABI£ II 


OXIDE PHASES PRESENT AFTER 


Alloy 

Time (hzrs) 

Hi Solid Soln. 

^^2°3 

Ni-13Cr-25Al-0.05Zr 

100 

VS 

m 


200 

vs 

s 


500 

vs 

w 


800 

vs 

w 

Ni-14Cr-24Al-0 . 07Zr 

800 

vs 

w 

Ni-14Cr-25Al-0. 12Y 

100 

vs 

m 


200 

vs 

w 


500 

vs 

m 


1000 

vs 

vs 

Ni-19Cr-24Al-0 . 04Zr 

100 

vs 

s 


500 

s 

vs 


800 

in . 

vs 


1600 

m 

vs 

Ni-17Cr-28Al-0 . 18Zr 

200 

s 

vs 


500 


s 


800 


m 


®ie oxide and metal phases were ranked by normalizing the largest peak 
sample. Ohe normalized peaks were ranked as 0. 8-1.0, very strong (vs) 
and 0 . 1 , very weak, (vw) . 


CTCLIC OXIDATION 


s 

pinel (a ) 

Other 

vw 

o 

(8.06-8.09) 


vw 

8.09 


in 

8.06-8.09 

CT 203 (vw) 

w 

8.06 



ro 

c» 


VW 

8.09 


vw 

s 

8.09-8.10 

(8.09) 

YA6 (vw) Yttrogamet 

(3Y203*5Al203) 

vw 

(8.06) 


w 

(8.06) 


vw 

(8.06) 


m 

(8.06) 


m 

(8,09) 

ZrO (w) monoclinic and cubic 
2 

vs 

(8.06) 

“2°3 Zr02(vw) 

vs 

(8.06) 

Cr^Os (w) 


for each phase to the largest peak measured on that particular 
0.5-0. 8, strong (s) j 0.2-0. 5, medium (m) } 0.1-0. 2, weak (w) ; 
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Figvire 6 Surface oxide morphologies after cyclic oxidation 

a. Ni-13Cr-25Al-0.05Zr alloy after 800 oxidation cycles 
(500x) b, Ni-14Cr-25Al-0. 12Y alloy after 1000 
oxidation cycles (lOOx) . 
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d 


Figure 6 cont. c. Ni-14Cr-24Al-0.07Zr alloy after 800 
oxidation cycles (lOOx) d. same as c, 
(500x) 
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e. Ni-19Cr-24Al-0. 04Zr alloy after 800 
oxidation cycles (200x) f. same as e. 
after 1200 oxidation cycles (lOOx) . 
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Figure 6 cont. g. Ni-19Cr-24Al-0.04 Zr alloy after 1600 oxidation 
cycles (lOOx) h. Ni-17Cr-28Al-0. 18Zr alloy after 
800 oxidation cycles (lOOx). 
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Figure 7 Surface oxide morphologies on the Ni-14Cr-25Al-0.12Y 

alloy after various cyclic oxidation exposures (lOOx) . 
a. 1 oxidation cycle b. 15 oxidation cycles 
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the micrographs are the locations where the oxide layer has spalled. After 
200 cycles, the oxide had spalled from over half the sample surface. No cor- 
relation. was found between the pattern of oxide spallation and the initial 

Y + 3 microstructure. After 1000 oxidation cycles (Figure 6b), the sample 
surface bears no resemblance to the early oxide scales. 

Oxidation of y + 3 NiCrAl alloys results in the growth of a near-surface 

Y layer by dissolution of the 6 phase at the y/y+ 3 interface. A1 is trans- 
ported across the y layer to the oxide/metal interface where it is consumed in 
formation of the oxide. To determine the width of the y layer, and to prepare 
for quantitative electron microscopy, the orthorhombic-shaped samples were 
sectioned perpendicular to the main axis (z or c), thereby exposing a cross 
section of the alloy perpendicular to the plane of the oxide scale. The cut 
samples were mounted and polished by standard metal lographic techniques. The 
near-surface Y-layei^ width was measured by optical microscopy using a filar 
eyepiece. The volume fraction of the 3 phase was also measured optically after 
various oxidation exposures. Following the optical examination, the mounted 
samples were carbon coated for use in the scanning electron microscope. 

Four to eight Ni , Cr, and A1 concentration/distance profiles were measured 
across the near-surface y layer for each alloy after various cyclic oxidation 
exposures. The concentration/distance profiles exhibited differences when 
measured across the y layer at various locations on a sample. The differences 
between the measured concentrati on/di stance profiles were usually greater than 
the possible discrepancies introduced by the measuring technique. The dif- 
ferences in the concentration profiles are therefore attributed to either com- 
positional variation within the as-cast sample, or to a varying requisite flux 
resulting from random oxide spallation. The concentration of the y and 3 
phases in the y + 3 region of each sample was also measured after various 



36 


oxidation exposures. Greater detail of the technique for measuring the con- 
centrations with the scanning electron microscope is given in Appendix A. The 
results will be presented separately for each alloy. 

Experimental Results 
Ni-13Cr-25Al-0.05Zr 

The near-surface y-phase layer thickened at a rapid, non-parabolic rate, 
the volume fraction of the 3 phase decreased from 0.73 to 0.43 after 800 
oxidation cycles. A third phase, y' > based on the Ni^Al-type intermetall ic 
phase formed at the y/y+3 interface between 500 and 800 oxidation cycles. In 
addition, the sample underwent a volume distortion after 500 oxidation cycles 
so that the sectioned plane (x-y plane shown in Figure 8) of the 800-hour sam- 
ple had decreased in the y direction but increased in the x direction (Figure 
8). Changes in the z direction could not be determined. The y layer, always 
measured in the x direction on all samples, was observed to increase abnor- 
mally due to the anomalous volume distortion. Considerable Kirkendal 1-type 
porosity was evident in the near-surface y layer (Figure 8a) but was not 
observed on other alloys exhibiting volume distortions (Ni-14Cr-25Al-0. 12Y 
alloy, Figure 8b). Therefore, the porosity alone cannot be the cause of the 
volume distortion. The microstructures and the near-surface y layer after 
various oxidation exposures of the Ni-13Cr-25Al-0.05Zr alloy are shown in 
Figure 9. 

The rapid loss of A1 results in a time-decreasing A1 concentration at the 
oxide/metal interface. The decreasing A1 concentration at this interface is 
attributed to an increased requisite flux encountered during cyclic oxidation 
accompanied by significant oxide spallation. Discontinuous loss of the oxide 
scale results in an average scale thickness less than that which occurs during 
parabolic oxidation. The thinner oxide scale results in a greater requisite 
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Original Width 


SOOpirt 


Figure 8 Volume distortion of sample coupons a. Ni-13Cr-25Al-0. 05Zr 
alloy after cyclic oxidation for 800 hours b. Ni-14Cr-25Al- 
0.12Y alloy after cyclic oxidation for 1000 hours. Schematic 
drawing to right shows the relationship between plane of sec- 
tioning (dashed) and sample coupon. 
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flux (JpgqOt l/^X^xide^' concentration at the oxide/metal interface 

decreases with time, increasing the concentration gradient, and thereby 
supplying the increased requisite A1 flux for continued AI 2 O 2 formation and 
growth. The decreasing A1 concentration at the oxide/nietal interface is shown 
as the left-most point of the A1 concentration/distance profiles in Figure 10. 
The Ni and Cr concentration profiles are also shown. The right-most point of 
each profile indicates the position of the y/y+ 3 interface, and the Ni , Cr, 
and A1 concentrations at that interface. The right-most point also gives the 
concentration of the y phase in the y+ 3 region of each sample. There were 
no measurable concentration gradients in the y+ 3 region. The error bars on 
the right-most point indicate the average standard deviation on the concentra- 
tions on that particular Ni , Cr, or A1 profile. Comparison of the individual 
profiles measured at various locations on a single sample indicated the scatter 
to be the result of profile-to-profile differences and not the result of the 
measuring technique. The differences between the concentration profiles 
measured at various locations along the sample are attributed to the non-uniform 
oxide spallation observed following cyclic oxidation of these alloys (Figure 7). 
The non-uniform oxide spallation results in a range in the oxide thickness on 
the sample surface, which further results in a range in the requisite flux 
along the sample surface. The differences in the measured concentration pro- 
files reflect the varying response of the alloy in supplying A1 over the range 
of requisite flux conditions. 

The A1 and Ni concentrations at the y/y + 3 interface (right-most concen- 
tration points. Figure 10) increase with increasing oxidation time, whereas the 
Cr concentration at the interface decreases with increasing oxidation time. In 
addition, the entire Ni concentration profiles shift to greater Ni concentra- 
tions with increasing oxidation time while the A1 and Cr concentration profiles 



Concentration (At.7o) 


0 



Distance (Mm) 


Figure 10. Ni , Cr, and A1 concentration/distance profiles in the 
Y layer after cyclic oxidation of the Ni-13Cr-25Al- 
O.OSZr and the Ni-14Cr-24Al-0,07Zr alloys. 
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generally decrease with increasing oxidation time. The overall increase or 
decrease in the concentration profiles can be observed and explained by 
plotting the concentration profiles as diffusion paths on the NiCrAl phase 
diagram. The diffusion paths are observed to "translate" toward the Ni-rich 
corner of the phase diagram, as shown in Figure 11. The solubility limit of 
A1 in the y phase (the line dividing the y and y+B fields) is the effective 
solubility limit at 1200°C used throughout this study, as discussed in 
Appendix B. 

The diffusion path translation can be understood by determining the com- 
position of the sample after oxidation. First, consider a Ni-18Cr alloy with 
the same Cr/Ni ratio (0.21) as the Ni-13Cr-25Al-0.05Zr alloy (point A in 
Figure 11). If pure A1 is added to this binary alloy, the resulting composi- 
tion of the NiCrAl alloy would lie on a line connecting point A to the 100 
at.%Al corner, also passing through the Ni-13Cr-25Al composition (see inset. 
Figure 11). Similarly, if A1 were removed from the Ni-13Cr-25Al alloy, the 
resulting composition would again fall on the dashed line in Figure 11. The 
composition of the Ni-13Cr-25Al -0.05Zr alloy after oxidation is the sum of the 
average compositions in the y+B region and the near-surface y layer, weighted 
in proportion to the fractional amount of the y+g region (or y layer) present 
after the oxidation exposure. The average composition of the y layer, deter- 
mined from Figure 10, and the average composition of the y+g region, deter- 
mined from the concentration and volume fraction of each phase, are indicated 
as open squares for each diffusion profile shown in Figure 11. The composition 
of the sample lies on a line connecting the average compositions of the y+g 
region and the y layer and is shown as open triangles in Figure 11. As shown, 
the composition of the sample after oxidation closely follows the dashed line 
of constant Cr/Ni ratio. 



100at.»/, 

Al 
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The diffusion path translation results as a natural consequence of losing 
A1 from the finite sample coupons during the oxidation process. Since the 
average. composition of the y layer resides to the left, or Cr-rich, side of 
the line of constant Cr/Ni ratio, the average composition of the y + 3 region 
must lie to the right, or Cr-poor, side of the line of constant Cr/Ni ratio. 

As the average composition of the y+g region shifts away from the Cr corner, 
the concentration in the y phase at the y/y+B interface also shifts in the 
same direction. The shifts in the concentration at the y/y + B interface tend 
to "pull" the diffusion paths in the y phase away from the Cr corner, resulting 
in the observed translation of the diffusion paths. The decrease in the 
volume fraction of the B phase with increasing oxidation time (Figure 9) is 
a result of the changing average composition in the y+B region associated 
with diffusion path translation. 

Diffusion path translation requires diffusion to occur in the y + B region 
of the alloy. It is apparent that the average composition of the y+B region 
increases in Ni concentration (Figure 10). Therefore, Ni must diffuse into 
the y+B region as A1 and/or Cr diffuse from the region. The fact that there 
were no measurable concentration gradients in the y+B region suggests that 
the region behaves as a thin, finite section; i.e., the average composition of 
the region changes uniformly with negligible concentration gradients. 

The Ni atoms must diffuse up the Ni concentration gradient in the y layer 
in order to enter the y+B region. A mass balance of Ni at the y/y+B inter- 
face, when considering the interface motion, requires the diffusion of Ni into 
the y+B region. The diffusion of Ni up its own concentration gradient is 
permissible when the Ni flux is significantly influenced by the concentration 
gra'^ients of either Cr or A1 . Diffusion of A1 down its ov/n concentration 
gradient in the y layer from the y+B region is supported by both the decrease 
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in the average A1 concentration in the y + 3 region (Figure 10) and by mass 
balance considerations regarding the motion of the y/y+3 interface. The 
interface motion also requires a flux of Cr atoms toward the moving interface 
(Appendix C). Since the average Cr concentration in the y + 3 region decreases 
with increasing oxidation time (Figure 11), part of the Cr supply comes from 
transport of Cr in the y+3 region to the y/y+3 interface. In addition, since 
the average Cr concentration in the y layer also decreases with increasing 
oxidation time (Figure 11), Cr must also diffuse down the Cr concentration 
gradient in the y layer to the y/y+3 interface. Therefore, changes in the 
average composition of both the y+ 3 region and the y layer, along with mass 
balance considerations at the y/y+3 moving interface, suggest that A1 and Cr 
diffuse down their respective concentration gradients while Ni diffuses up its 
own concentration gradient. 

Diffusion of Ni away from the oxide/y-layer interface is not contradicted 
by the observation of mixed Ni oxides in the scale (Table II). Considerable 
oxide cracking or spalling to the metal substrate would allow a high concen- 
tration of oxygen to contact the metal. Significant amounts of NiO may 
initially form and rapidly grow since Ni is the most abundant species at the 
oxide/metal interface. As A1 diffuses to the interface, the A1 will react 
with the NiO to form the more thermodynamically stable NiAl20^. When the 
NiO is completely converted, an Al^O^ scale will form below the NiAl^O^. 

Figure 12 shows a duplex scale structure for the Ni-13Cr-25Al-0.05Zr alloy 
after 800 hours of cyclic oxidation. The area indicated A is the metal sub- 
strate containing obvious oxide grain imprints. The rough area designated B 
was determined by energy-dispersive spectroscopy (EDS) to contain roughly 60% 
Cr, 30% Ni , and 10% A1 , probably a mixture of NiCr^O^ and NiAl20^ spinels. 

The area designated C shows columnar grains containing mainly A1 , most 



Figure 12. Scanning electron micrograph of the surface oxide scale on 
the Ni-13Cr-25Al-0.05Zr alloy after 800 oxidation cycles, 
a. oxide scale (dark) and metal substrate where oxide has 
spalled (light) , (600x) ; b. high magnification of oxide 
scale edge showing duplex oxide structure; note oxide im- 
prints in the metal substrate, A; mixed spinel oxides, B; 
and columnar Al^O^ grains, C. (2400x) 
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probably Al^O^ (oxygen could not be detected with the EDS), growing below the 
spinels into the rnetal substrate. The explanation of initial NiO formation 
converting to aspinel with subscale formation of Al^O^ is similar to models 
of transient oxide formation during isothermal oxidation (5,21,23). 

Ni-14Cr-24Al-0.07Zr 

The oxide scale provided by this alloy affords greater oxide protection 
during cyclic oxidation than the previously discussed Ni-13Cr-25Al-0.05Zr alloy 
with approximately the same composition. The present O.OZZr alloy was alloyed 
with Zr by induction melting elemental Ni , Cr, and A1 in zirconia crucibles, 
allowing the Zr to be leached from the crucible. Two points are immediately 
obvious in a comparison of the microstructures and concentration profiles of 
each alloy. First, the y layer of the O.OZZr alloy is much less than the pre- 
viously discussed 0.05Zr alloy, both after 800 oxidation cycles (Figure 13). 
Secondly, the A1 concentration at the oxide/metal interface after 800 oxidation 
cycles is considerably greater for the O.OZZr alloy than that for the O.OSZr 
alloy after 800 or even *500 oxidation cycles. The concentration/distance 
profiles for the O.OZZr alloy are shown as open points in Figure 10. The 
concentration profiles (Figure 10), taken together with the difference in the 
y-layer widths (Figure 13) and the weight change data (Figure 4), indicate that 
the O.OZZr alloy possesses much better spalling characteristics than the O.OBZr 
alloy. A study of the effect of the Zr concentration on the cyclic oxidation 
of NiCrAlZr alloys by Barrett, Kahn, and Lowell (48) indicates that the dif- 
ference in the Zr content of the two alloys does not explain the variation 
observed in 'their oxidation behavior. Hence it appears that the variation in 
the cyclic oxidation behavior may be attributable to the method of alloying 
the Zr and that Zr pick-up from the zirconia crucibles during the induction 
melting process improves the cyclic oxidation behavior of these alloys. The 
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Isothenral 



1000 hrs. Ni-14Cr-24Al-0.07Zr 


Cyclic and Isothermal oxidation 
of Ni-13Cr-25Al-0.05Zr alloy and 
Ni-14Cr-24Al-0.07Zr alloy. All 
micrographs at lOOx. 
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improved cyclic oxidation behavior resulting from Zr pick-up from the zirconia 
crucibles is under further investigation at the NASA-Lev/is Research Center (57), 

The Ni-14Cr-24Al-0.07Zr alloy was also examined after isothermal oxidation 
for 200 and 1000 hours at 1200°C. The y-layer width after isothermal oxidation 
is very slight in comparison to the width following cyclic oxidation. The y- 
layer widths and resulting microstructures after isothermal and cyclic oxidation 
of the 0.07Zr alloy and the cyclic oxidation of the 0.05Zr alloy are compared in 
Figure 13. The 3 phase in both alloys underwent considerable coarsening during 
the high-temperature exposure. The volume fraction of the 3 phase remained 
essentially constant during isothermal oxidation (0.64 - 0.68), whereas the 
volume fraction of the 3 phase significantly decreased during cyclic oxidation. 
The difference in the volume fractions of the 3 phase is directly attributable 
to the lack of any significant diffusion path translation in the isothermal ly 
oxidized alloy. The relatively large loss of A1 during cyclic oxidation re- 
sults in significant diffusion path translation, changes in the average compo- 
sition of the y+3 region, and consequently, a change in the volume fraction 
of the 3 phase. 

Limited concentration measurements in the y layer after isothermal' oxida- 
tion up to 1000 hours indicate a very low A1 concentration gradient. A low A1 
concentration gradient is consistent with the expectation of a low requisite 
flux during parabolic oxidation. The near-surface y layer was measured to 
thicken parabolically with time. The thin black oxide stringers growing in 
from the oxide scale are common in Zr-containing alloys (27, ,48). The growth 
of the oxide stringers into the substrate appears irrespective of the presence 
of either the y or 3 phase (Figure 14). EDS analysis of the stringers indicates 
a high A1 content (presumably Al^O^) with a very slight trace of Zr. 




Figure 14. Growth of the ^1^0 stringers into the 
substrate of the Ni-14Cr-24Al-0 . 07Zr 
alloy after 1000 hours of isothermal 
oxidation. (300x) 
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Ni-19Cr-24Al-0.04Zr 

This alloy exhibited the best cyclic oxidation resistance of the four 
alloys studied. Strong Al^O^ formation (Table II) resulted in a low loss of 
scale during cyclic oxidation, as indicated by the weight change data in Figure 
4. The near-surface y layer increased at a much lower rate than the other 
cyclically oxidized alloys. Figure 15 shows the alloy microstructures after 
cyclic oxidation to 1600 hours. The volume fraction of the 3 phase decreased 
only slightly from approximately 0,62 to 0.53 after 1600 oxidation cycles. 

The good cyclic oxidation behavior of this alloy is reflected in the low 
A1 concentration gradients in the y layer. A low requisite A1 flux associated 
with good oxide scale adherence results in a very small decrease in the A1 
concentration at the oxide/metal interface. The concentration/distance profiles 
are shown in Figure 16. The 800-hour A1 concentration profile appears abnor- 
mally high, most probably the result of an increased A1 content produced at 
the casting stage. The relatively low loss of A1 from this alloy results in a 
small amount of di ffusion path translation, shown in Figure 17. The difference 
in the concentration profiles resulting from a low, and high, requisite flux 
condition can be seen by comparing the concentration profiles of the 
Ni-13Cr-25Al-0,05Zr alloy and the Ni -19Cr-24Al -0 .04Zr alloy, shown in Figures 
10 and 16, respectively. The obvious di f ference between the concentration 
profiles is the A1 concentration at the oxide/metal interface. The interfacial 
A1 concentration of the Ni - 19Cr-24Al -0 . 04Zr alloy remains relatively high 
(> 8 at.% after 1600 hours) , whereas the interfacial concentration of the 
Ni-13Cr-25Al-0.04Zr alloy is considerably lower (< 7.5 at.% after only 800 
hou rs . 

The good oxide adherence of Ni-19Cr-24Al -0.04Zr al loy may be due to 
accommodation of the thermal ly induced strains by the increasingly convoluted 
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Figure 15. Micrographs of the near-surface Y layer after cyclic 
oxidation of the Ni-19Cr-24Al-0.04Zr alloy. 
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oxide/metal interface. Oxide ridges are apparent in Figures 6c- 6g. The con- 
voluted oxide/tnetal interface is also apparent in the micrographs of Figure 15 
(800- 1600 hour exposures). The source containing References 13 and 15 dis- 
cusses the theoretical and experimental aspects of stress effects and oxide 
morphologies during oxidation. 

Ni-14Cr-25Al-0.12Y 

The Y content of this alloy effects cyclic oxidation behavior very simi- 
lar to that of the Ni -13Cr-25Al -0.05Zr alloy. The near-surface y layer 
increases at approximately the same rate for both alloys through 500 hours. 

The difference at longer times is probably due to different volume distortions 
which occur in both alloys at the longest oxidation time (Figure 8b). The 
similarity in the y-layer width between the two alloys can be observed by 
comparing the micrographs of the Y-containing alloy (Figure 18) with those of 
the Zr-containing alloy (Figure 9). 

The concentration/distance profiles also exhibit characteristics similar 
to the Zr-containing alloy. The A1 concentration at the oxide/metal interface 
significantly decreases with increasing oxidation time, while the A1 concen- 
tration at the y/y+3 interface similarly increases as a result of diffusion 
path translation. The concentration/distance profiles are shown in Figure 19, 
Two A1 concentration profiles are shown for the 1000-hour sample near the 
oxide/metal interface. The lower A1 concentration profile resulted from 
measuring the concentrations in the metallic sections which extended past' the 
average oxide/metal interface position. The diffusion path translation is 
shown in Figure 17. 

Two important points are worth noting in a comparison of the Y- and Zr- 
containing alloys. First, the A1 concentration at the y/y+B interface is 
approximately the same for both alloys at each oxidation time. This similarity 
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implies that Y does not significantly affect the maximum solubility of A1 in 
the Y phase. It is also shown in Appendix B that Zr does not appear to signi- 
ficantly affect this solubility. Second, although both alloys have approxi- 
mately the sane y-layer width after equivalent oxidation times, the A1 concen- 
tration gradient in the Y-containing alloy is greater than that in the Zr- 
containing alloy. The similar y-layer widths imply that both alloys have lost 
A1 at approximately the same rate, further implying that both alloys have 
approximately the same requisite flux. The difference in the concentration 
gradients, but the same requisite flux, suggests that Y adversely affects the 
supply of A1 to the oxide/metal interface. Since Appendix C suggests that the 
A1 flux in the y layer is mainly dependent on the A1 concentration gradient, it 
may be inferred that Y decreases the A1 diffusivity below that in the Zr- 
containing alloy. A greater concentration gradient is therefore required in 
the y layer to supply the same requisite flux as the Zr-containing alloy. 

The oxide stringers observed after oxidation of Y-containing alloys are 
generally attributed to the selective oxidation of grain-boundary yttrides. 

Figure 20 shows the selective oxidation of a grain boundary yttride after lOQ 
hours of cyclic oxidation of the Ni-14Cr-25Al-0. 12Y alloy. EDS analysis of the 
unoxidized portion of the yttride indicated a composition of mainly Ni with 
lesser amounts of Y and trace amounts of A1 and Cr. 

Ni-17Cr-28Al-0.18Zr 

The high Zr content of this alloy effects poor oxide spalling character- 
istics similar to the Ni -13Cr-25Al -0.05Zr and Ni -14Cr-25Al -0. 12Y alloys. The 
A1 content, being the highest in this alloy, results in a lower growth rate 
of the near-surface y layer while still supplying A1 to the oxide scale at 
approximately the sarre rate as the Ni-13Cr-25Al-0.05Zr and the Ni-14Cr-25Al -0. 12Y 
alloys. The width of the y layer in the Ni-17Cr-28Al-0. 18Zr alloy is shown in 
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Figure 20. Selective oxidation of grain boundary yttrides 
on the Ni-14Cr-25Al-0. 12Y alloy after cyclic 
oxidation for 100 hours. (50Qx) 
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Figure 21. The high Zr content also results in an increased occurrence of 
Al^Og stringers in the y layer, apparent in Figure 21. The increased number 
of oxide stringers accounts for the abnormally large, positive weight change 
during cyclic oxidation (Figure 4). 

More interesting in the Ni-17Cr-28Al-0. 18Zr alloy is the occurrence of 
localized breakaway oxidation after 500 and 800 oxidation cycles. Internal 
oxidation of A1 , signaling external formation of less protective oxides of Ni 
and Cr, is shown in Figure 22. An increased slope in the weight change data 
after approximately 500 oxidation cycles (Figure 4), indicating rapid spalla- 
tion of poorly adherent scales, further attests to the onset of breakaway 
oxidation. Breakaway oxidation occurs in localized regions because of the non- 
uniform nature of oxide spallation from the sample surface. The x-ray maps of 
Figure 23 show Cr-rich oxides (probably Cr^O^) forming in an internally oxi- 
dized subscale region indicating significant oxygen influx into the alloy. 
Figure 22 also shows a partial link-up of Al^O^ precipitates at the boundary 
between the y layer and the internally oxidized zone. 

Breakaway oxidation results from a failure of the alloy to supply the 
requisite A1 flux for sole formation and growth of an Al^O^ oxide scale. As 
was evident from the previously discussed alloys, the A1 concentration at the 
oxide/metal interface decreased with time, increasing the A1 concentration 
gradient in the y layer, and allowing the alloy to supply the requisite flux 
to the oxide scale. Decrease of the interfacial A1 concentration to zero 
prevents further supply of the requisite flux, and less-protective oxides will 
eventually form. Formation of less-protective oxides permit a high influx of 
oxygen into the alloy substrate internally oxidizing the A1 . Further oxida- 
tion results in either total internal oxidation of the A1 in the sample, with 
rapid spalling of Ni and Cr oxides, or internal oxide link-up resulting in a 
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Figure 22. Localized breakaway oxidation of the Ni-17Cr-28Al--0. 18Zr 
alloy after 800 oxidation cycles a. (30x) b. Internal 
oxidation of A1 and some Cr. (180x) 




63 



Figure 23, X-ray maps of internal oxidation substrate shown in 
Figtxre 22. a. Internally oxidised subscale (180x) 
b. A1 x-ray map. 
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new AI^O^ protective barrier. The greater the volume fraction of internally- 
oxidized precipitates deposited, the more tortuous the path the incoming oxygen 
must travel and the greater the chance of precipitate link-up establishing a 
continuous, protective oxide scale. Rapp (58) has reviewed the theoretical 
and experimental work regarding the solute concentration and diffusivity 
necessary for the transition from internal oxidation to the formation of an 
external, continuous oxide scale. It is apparent from Figure 23 that the 
supply of A1 to the internal oxidation/y-layer interface is sufficient to 
reestablish a continuous Al^O^ barrier by link-up of the Al^O^ precipitates. 

The A1 concentration at the oxide/metal interface for the Ni-17Cr-28Al- 
O.lSZr alloy decreases to nearly zero after 800 oxidation cycles. The A1 
concentration profiles showed a significant difference in concentrations 
between the five to eight profiles measured at various locations along the 
sample surface. The profiles are shown in Figure 24. The two error bars 
below the A1 concentration profiles are indicative of the differences between 
various concentration profiles measured on the same sample. The concentrations 
at the y/y+ 3 interface change in accordance with the diffusion path transla- 
tion shown in Figure 17. 

The occurrence of breakaway oxidation on the 500 and 800 hour samples 
results from a failure of the alloy to supply the requisite flux when the A1 
concentration at the oxide/metal interface decreased to nearly zero. It is 
interesting to note that this alloy with the highest A1 content (28 at.%) 
undergoes breakaway oxidation before the alloys with less A1 , even though the 
requisite A1 flux and the weight of A1 consumed in this alloy will be shown to 
be similar to that in the Ni-13Cr-25Al-0.05Zr and Ni-14Cr-25Al -0. 12Y alloys. 
This point will be discussed further in the next section. 
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. Ni, Cr, and A1 concentration/distance profiles in the 
Y layer after cyclic oxidation for the Ni-17Cr-28Al- 
O.lSZr alloy. 


Figure 24 
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Summary 

The observations concerning the A1 transport during the cyclic oxidation 
of NiCrAlZr(Y) alloys will be briefly reviewed. The differences between iso- 
thermal and cyclic oxidation will also be noted. The observations from the 
experimental analysis are: 

1. Oxide spallation during thermal cycling results in a more rapid loss 
of A1 over that which occurs during isothermal oxidation. The accel- 
erated loss of A1 is due to the higher requisite flux associated with 
cyclic oxidation. 

2. Repeated oxide spallation results in a decreasing A1 concentration at 
the oxide/y-layer interface. The interfacial A1 concentration de- 
creases, increasing the concentration gradient, in order to supply 
the higher requisite flux. No measurable decrease in the interfacial 
A1 concentration was observed during isothermal oxidation. 

3. The rate of A1 decrease at the oxide/y-layer interface was in propor- 
tion to the severity of the oxide spallation. The Ni-13Cr-25Al-0.05Zr, 
the Ni-14Cr-25Al-0.12Y, and the Ni-17Cr-28Al -0. 18Zr alloys, all with 
poor spalling characteristics, exhibited large decreases in the inter- 
facial A1 concentration after 800 - 1000 hours. In comparison, the 
Ni-19Cr-24Al-0.04Zr alloy, with good spalling characteristics, exhifi- 
ited a small decrease in the interfacial concentration after 1600 
hours. 

4. Breakaway oxidation was observed when the A1 concentration at the 
oxide/y-layer interface decreased to approximately zero in the 
Ni-17Cr-28Al-0. 18Zr alloy. Breakaway oxidation resulted from a 
failure of the alloy to continue to supply the requisite flux. 

The y-layer width following cyclic oxidation was considerably greater 


5, 
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than that observed after isothermal oxidation. The growth rate of 
the Y layer in the cyclically oxidized alloys reflects the influence 
of the increased requisite flux associated with oxide spallation. 

6. The loss of A1 from finite samples led to diffusion path translation, 
changes in the composition of both phases in the y+B region, and a 
decrease in the volume fraction of the 3 phase. The isothermally 
oxidized alloys, having lost relatively little A1 , showed insignifi- 
cant diffusion path translation or decrease in the volume fraction 
of the 3 phase. 
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NUMERICAL MODELING RESULTS AND DISCUSSION 


Numerical Modeling 

Numerical models may be used to simulate and predict the diffusional 
transport in alloys at elevated temperatures. Most numerical models are based 
on finite-difference techniques which permit considerable flexibility when a 
closed-form analytical solution is not possible. The flexibility of a numeri- 
cal model allows incorporation of concentration-dependent diffusion coefficients, 
complex boundary conditions, non-isothermal temperature conditions, and over- 
lapping diffusion fields. A number of binary and ternary numerical models 
utilizing finite-difference techniques have been developed to simulate diffu- 
sion under various conditions (59-62). Even with the increased flexibility, 
numerical models still require some assumptions in order to determine a 
solution. 

Numerical modeling of atomic diffusion in systems containing moving 
boundaries requires solution of Pick's second law in all single-phase regions 
constrained by mass balances at the moving boundaries. If concentration- 
independent diffusion coefficients are assumed. Pick's second law for the 
NiCrAl system can be stated as (63): 


dC 


A1 

dt 


"AlAl 


.2 ^AlCr 


dx‘ 


dx‘ 


(la) 


(lb) 


where C^-j and are the A1 and Cr concentrations, respectively. The parameters 
X and t refer to distance and time, respectively. The four diffusion coeffi- 
cients— D^^^p and — 


will be defined shortly. 
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The mass balance at the y/y + 3 interface has the form: 
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where: ^ = velocity of the y/y+ 6 interface. 


A1 flux in the y phase evaluated at the y/y+B interface. 


jY 

^A1,5 

= A1 concentration in y+B region. 


(2a) 


(2b) 


^Y - 
■'Al 


A1 concentration in the y phase at the y/y+B interface. 


The Cr concentration parameters in Eq. 2b are defined similarly to those for 
the A1 . The A1 and Cr fluxes are defined by the right-most term within the 
brackets in Eqs. 2a and 2b. afid *^crCr ^ ^ 

direct diffusion coefficients relating the A1 concentration gradient to the A1 
flux and the Cr concentration gradient to the Cr flux (right-most terms, Eqs. 
2a and 2b), respectively. and are the cross term, or indirect, 

diffusion coefficients relating the Cr concentration gradient to the A1 flux 
and the A1 concentration gradient to the Cr flux (right-most terms, Eqs. 2a 
and 2b), respectively. The terms containing and in Eqs. 1 and 2 

are referred to as the cross terms which describe the effect of one component 
on the diffusion of other components in ternary and more complex systems (63). 

Loss of a component from the metal substrate at a free surface results in 
recession of the metal substrate. The appropriate equation describing the 
metal recession due to loss of A1 from the y phase in NiCrAl alloys is: 


^A1 


dC 


2 _ 


A1 dt 
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(3a) 
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where: ^ = velocity of the metal substrate recession. 


'A1 


partial molar volume of A1 in the y phase. 


C ^1 = A1 concentration in the y phase at the oxide/y- layer interface, 

= A1 flux in the y phase evaluated at the oxide/y-layer interface 
(defined similarly to the right-bracketed term in Eq. 2a). 

The greater the partial molar volume of A1 , the greater the surface recession 
caused by loss of A1 from the metal substrate. 

The second goal of this research consisted in developing a model to simu- 
late the A1 diffusion in y+3 NiCrAl alloys during cyclic oxidation. As pre- 
viously shown, the diffusion of A1 to form an AI 2 O 2 oxide scale results in the 
growth of a near-surface y layer. The loss of A1 from the metal substrate 
resulted in the recession of the Al^O^/y-layer interface. Hence, it was 
necessary to determine the A1 diffusion in the near-surface y layer constrained 
by the y/y + B and Al202/y-layer moving boundaries. The cross term effect of 
the Cr concentration gradient was neglected on the basis that was much 

greater than °A1 as suggested by the diffusion couple work in Appendix C. 
Therefore, the A1 transport was modeled on the assumption that the A1 concen- 
tration gradient was the dominant driving force for A1 diffusion. 

A numerical model was developed capable of predicting A1 concentration/ 
distance profiles, the y/y+B interface motion, the oxide/y-layer surface 
recession, and the weight of A1 consumed. In addition, the model could pre- 
dict the time to breakaway oxidation by predicting the time the A1 concentra- 
tion at the oxide/y-layer interface approached zero. Seven assumptions were 
made in the development of the numerical model. They were: 

1. The diffusion coefficients were independent of concentration. 

2. The y/y+B interface was planar. 

3. The oxide/y-layer interface was planar. 
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4. The partial molar volume of A1 was independent of concentration and 
phase. 

5. Sample heat-up and cool-down times were rapid; i.e., diffusion only 
occurred at 1200°C. 

6. Only Al^O^ was formed. 

7. The driving force for A1 was solely the result of the A1 concentration 
gradient. 

The finite-difference forms of Pick's second law and the mass balances at 
each interface were utilized to describe the A1 diffusion. The resulting 
diffusion equations to be solved were: 




(4a) 


in the y layer; 




- C°) ^ = J] 


= -D 


dx 


at the y/Y + B interface; and 


- ''ai) .s ^^2 ^ ,Y 


V 


A1 


dt ^2 


= -D 


dx 


(4b) 


(4c) 


at the oxide/y-layer interface, where the concentration and flux parameters 
refer only to A1 , and D refers to the diffusion coefficient approximating 
as discussed in Appendix C. All other parameters are described previously for 
Eqs. 1, 2, and 3. 

The boundary condition at the oxide/Y-layer interface may be stated as 
either a time-dependent flux (Jp or as a time-dependent interface concentra- 
tion (C ). If the flux (J^) or the specific rate of A1 consumption (i.e., the 
requisite flux) is known, the interfacial A1 concentration can be determined 
from the right-most term in Eq. 4c and Pick's second law in the y layer (Eq. 4a). 
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Similarly, if the A1 concentration at the interface (C^) is known, Eqs. 4a and 
4c alloy the flux (Jp to be determined. The interfacial A1 concentration (C^), 
estimated from the concentration/distance profiles measured on the cyclically 
oxidized alloys, was input to the model. The range of concentration gradients 
measured on the cyclically oxidized alloys did not facilitate easy estimation 
of the flux (Jp. 

The time dependence of the A1 concentration at the y/y+ 3 interface was 
also input to the model. The A1 concentration at this interface increased with 
time as the diffusion paths translated toward the Ni-rich portion of the phase 
diagram. The time dependence of the A1 concentration at the y/y + 3 interface 
and at the oxide/Y-layer interface is shown in Figure 25. The finite-difference 
form of the diffusion equations and the basic operations of the numerical model 
are described in greater detail in Appendix D. The determination of the dif- 
fusion coefficients for each alloy is shown and discussed in Appendix C. 

Results and Discussion 

A1 Concentration Profiles and y-Layer Widths 

The numerical model is able to predict accurately the A1 concentration 
profiles and y-layer widths through 500 hours for the Ni -13Cr-25Al -0. 05Zr alloy. 
The anomalous volume expansion of this alloy results in the underestimation of 
the Y-1ayer width for the 800-hour prediction. The predicted and measured A1 
concentration/distance profiles are shown in Figure 26. 

Similarly, the concentration/distance profiles and y-layer widths are 
accurately predicted for the Ni-14Cr-25Al-0. 12Y alloy through 500 hours. Again 
the anomalous volume expansion results in the underestimation of the Y“^aye>" 
width at 1000 hours. The concentration/distance profiles are shown in Figure 
27. The curvature in the 500-hour and 1000-hour A1 concentration profiles is 
greater than the curvature in the predicted profiles. The additional curvature 
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Figure 26. Predicted and measured Al concentration/distance profiles for the 
Ni-13Cr-25Al-0.05Zr alloy. 
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Figure 27. Predicted and measured Al concentration/distance profiles for the 
Ni-14Cr-25Al-0.12Y alloy. 
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may be the result of either a concentration dependence of the diffusion coeffi- 
cient, the diffusion coefficient decreasing as the A1 concentration decreases,^ 
or the effect of the Cr concentration profile on the A1 diffusion. 

The concentration/distance profiles and y-iayer widths were predicted 
with moderate success for the Ni-19Cr-24Al -0.04Zr alloy. This alloy exhibited 
the best oxide spalling characteristics of the four alloys studied and, conse- 
quently, has lost the least A1 to the formation of Al^O^- Since the requisite 
flux for this alloy was relatively low, the concentration gradients measured 
and predicted were also low in comparison with the other alloys. The measured 
and predicted concentration/distance profiles for this alloy are shown in 
Figure 28. Small changes in the time dependence of the A1 concentration at 
the oxide/metal interface, or at the y/y+ 3 interface, resulted in large changes 
in the y-layer width. In addition, small differences in the initial alloy com- 
position could also result in large differences between the predicted and 
observed behavior. An example of this is the high A1 concentration measured 
in the y layer of the 800-hour sample. The predicted concentration profile is 
much lower than that measured resulting in an overestimation of the y-layer 
wi dth . 

The concentration/distance profiles and y-layer width for the Ni-17Cr-28Al- 
0.18Zr alloy were also predicted with moderate success. The width of the y 
layer was predicted fairly well through 500 hours but overestimated at 800 
hours. The 500-hour and 800-hour samples were both from a different melt than 
the shorter-time samples and appear to have a slightly higher A1 concentration 
at the y/y+3 interface than what would be expected based on the shorter-time 
concentration profiles. The measured and predicted concentration/distance 
profiles are shown in Figure 29. 

The longer-tine concentration profiles exhibit more curvature than the 
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predicted concentration profiles. The additional curvature could be due to 
either concentration-dependent diffusion coefficients or the cross-term effect 
of the Cr concentration profiles. The curvature in the concentration gradients 
near the oxide/metal interface may also be due to the very irregular shape of 
this interface caused by the protrusion of the AI 2 O 2 stringers into the y layer. 
The A1 diffusion in this region may not be in the single direction tov/ards the 
oxide/metal interface but in three directions as A1 also diffuses to the oxide 
stringers in a direction parallel to the sample surface. The increased sinks 
for the A1 diffusion could effect lower A1 concentrations and steeper A1 
gradients in the regions containing the oxide stringers. 

The ability of the model to predict the width of the y layer for each of 
the four cyclically oxidized alloys is shown in Figure 30. The Ni -14Cr-24Al - 
0.07Zr alloy after isothermal oxidation is shown also. Good agreement is 
evident for each of the alloys excluding the 800-hour Ni-13Cr-25Al-0.05Zr sam- 
ple and the 1000-hour Ni-14Cr-25Al-0.12Y sample, both of which exhibited 
anomalous volume expansions. 

The agreement between the predicted and experimental ly determined weight 
of A1 consumed is shown in Figure 31. The experimental values were determined 
graphically by the use of Figures 26 - 29. The anomalous volume expansions 
and the inability to determine accurately the oxide/metal surface recession at 
longer times prohibited the experimental determination of the weight of A1 con- 
sumed beyond 500 hours. Because of the scatter in the concentration values in 
Figures 26 - 29, the accuracy of the experimental values is estimated to be 
± 15%. Figure 31 shows that the Ni-19Cr-24Al -0. 04Zr alloy exhibits the lowest 
loss of A1 , consistent with the observation of good spalling characteristics 
attributed to an adherent AI 2 Q 2 oxide scale. The alloys Ni-13Cr-25Al-0.05Zr, 
Ni-14Cr-25Al-0. 12Y, and Ni-17Cr-28Al-0. 18Zr exhibit a high loss of A1 indicating 
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Figure 31. Predicted and measured weight of A1 consumed. 
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poor scale adherence during cyclic oxidation. It is interesting to note that 
the alloy with the best scale adherence contained the lowest Zr content which 
was obtained by leaching Zr from a zirconia crucible during the melting process. 

Breakaway Oxidation 

Predicting the time to breakaway oxidation can be illustrated for the 
Ni-13Cr-25Al-0.05Zr alloy and the Ni- 17Cr-28Al -0. 18Zr alloy. Figure 31 indi- 
cates that both alloys have lost approximately the same amount of A1 after 800 
hours of cyclic oxidation. Hence, the requisite fluxes for both alloys were 
approximately equal. After several hundred oxidation cycles, the requisite 
flux approaches a constant value as the average oxide scale tends toward a 
constant thickness. The constant oxide-scale thickness is a consequence of 
scale growth being equal to scale spallation during each cycle. This behavior 
results in the linear portion of the paral inear weight change curves shown in 
Figure 4. 

An approximately equal time-independent flux (Jp was input to the numeri- 
cal model as a boundary condition at the oxide/y-layer interface for the 
Ni-13Cr-25Al-0.05Zr and Ni-17Cr-24Al-0. 18Zr alloys. The A1 concentration at 
the oxide/y-l ayer interface of the higher-Al -containing alloy was predicted to 
approach zero at a much shorter time. The input flux and the predicted inter- 
facial A1 concentration is shown in Figure 32. A parabol ical ly decreasing flux 
and the resulting time-independent interfacial A1 concentration for the iso- 
thermal oxidation of the Ni-14Cr-24Al -0.07Zr alloy is also shown. The experi- 
mental scatter in the A1 concentration near the oxide/y-layer interface in the 
Ni-17Cr-28Al-0. 18Zr alloy suggests that breakaway oxidation over much of the 
sample surface could occur before 1000 oxidation cycles. The observation of 
limited breakaway oxidation after 500 hours, and the occurrence of greater 
areas undergoing breakaway oxidation after 800 hours (Figure 22), confirm the 
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Figure 32. 


Input flux (J _ ) and the resulting A1 concentration 

surface 

at the oxide/metal interface. The flux and interfacial 
A1 concentration for isothermal oxidation is also shown. 
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model's prediction of the early failure of the high-Al alloy even though both 
alloys have lost approximately equal amounts of Al. 

The time to breakaway oxidation is dependent on the alloy composition and 
the Al diffusivity. The Al content, as just shown, is not a sufficient param- 
eter for ranking alloys according to their ability to resist breakaway oxida- 
tion. The alloy composition, specifically with regard to its position in the 
y+3 region of the phase diagram, and the diffusion coefficient, relating the 
Al flux to the concentration gradient, must also be considered when predicting 
breakaway oxidation. The importance of the diffusivity and alloy position in 
the phase diagram can be seen in a comparison of the two alloys previously 
discussed, viz., the Ni-13Cr-25Al-0.05Zr and the Ni -17Cr-28Al -0. 18Zr alloys. 

The diffusivity of the low-Al Ni-13Cr-25Al-0.05Zr alloy was measured as twice 
that of the high-Al Ni-17Cr-28Al-0. 18Zr alloy. Therefore, the low-Al alloy 
could supply the same Al flux to the oxide/metal interface as the high-Al alloy 
but with only half the Al concentration gradient. The difference in the con- 
centration gradients near the oxide/metal interface is obvious from a compari- 
son of the predicted Al concentration profiles for each alloy at 800 hours, 
shown in Figure 33, where both alloys are supplying the same requisite flux 
(Figure 32). 

Two important points regarding the bulk alloy composition and position in 
the phase diagram can be noted in Figure 33. First, the high-Al alloy has a 
large concentration drop at the y/y+ 3 interface, from the bulk concentration 
(Cq) to the maximum solubility of Al in the y phase (C^g)- The large concen- 
tration difference at the y/y+3 interface allows the high-Al (28 at.%Al) alloy 
to supply the same amount of Al to the oxide after 800 hours with less than 
three-fourths the depletion width of the low-Al (25 at.%Al) alloy (see Figures 
lb and Ic). The second point to note is that the position of the high-Al alloy 
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Al 




Figure 33. Predicted Al concentration profiles for 800 
oxidation cycles of the Ni-17Cr-28Al-0. 18Zr 
and the Ni-13Cr-25Al-0.05Zr alloys. 
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in the y + 3 field results in a low value of the maximum solubility of A1 in 
the Y phase (C^g)* Consequently, the A1 concentration at the oxide/y-layer 
interface is lower than that in the low-Al alloy. The same statement would be 
true even if the diffusivities, and therefore the concentration gradients, of 
both alloys were equal (see Figures lb and Id). An optimum y + 3 alloy compo- 
sition may be selected to maximize the time before breakaway oxidation. The 
alloy would be chosen with the highest possible A1 solubility in the y phase 
(C^g) and the maximum A1 content (C^) effecting a low y-layer growth rate. An 
alloy fitting this description is indicated as point A on the ternary phase 
diagram in Figure 33. 

The assumptions made in the development of the numerical model can be 
examined to identify some of the possible sources of inaccuracy introduced to 
the model. The assumption of a planar interface at the y/y + 3 boundary appears 
justified by the observation that the boundary remains planar on a macroscopic 
scale (Figures 9, 13, 15, 18, and 21) and by the fact that no concentration 
gradients were observed in the y + 6 region of any sample. The assumption of a 
constant molar volume seems reasonable based on the observation that consider- 
able variation in the partial molar volume results in relatively small changes 
in the overall diffusion behavior (Appendix D). The assumption of exclusive 
formation of Al^O^ was accepted even though some oxides containing Ni and Cr 
were observed. Formation of Ni or Cr oxides would accelerate the motion of the 
oxide/metal interface but v'oulo be expected to have no greater effect on the 
overall diffusion behavior than the partial molar volume. The assumption that 
negligible diffusion occurs during heating and cooling of the sample appears 
justified considering the exponential relationship of the diffusion coefficient 
with temperature. The assumption of a planar interface at the oxide/y-layer 
interface was only an approximation to the actual case, especially for the alloy 
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containing the highest Zr concentration (Ni-17Cr-28Al-0. 18Zr) . This assumption 
may have accounted for some of the inacuracy in the predicted concentration 
profiles especially near the oxide/y-layer interface. The assumption of negli- 
gible cross-term effects (i.e., >> ^AlCr^ adequately discussed 

in Appendix C. No conclusive justification for the concentration independence 
of the diffusion coefficients can be given. Any concentration dependence of 
the diffusion coefficient could easily account for most of the discrepancy 
between the predicted and measured concentration profiles. Based on the dif- 
fusion data for the similar CoCrAl systems (64), could easily vary with 

composition by an order of magnitude. A further assumption, not explicitly 
stated, was that the initial composition of each sample, oxidized at different 
times, was the same. Possible exceptions to this assumption were noted for 
alloys Ni-19Cr-24Al-0.04Zr and Ni-17Cr-28Al -0. 18Zr. 


Summary 

The results of numen'cally modeling the A1 transport during the cyclic 
oxidation of y+3 NiCrAl alloys will be briefly reviewed. The results were: 

1. The A1 concentration/distance profiles resulting from cyclic oxida- 
tion of the y+B alloys were predicted fairly accurately through 500 
hours. The anomalous volume distortions in two of the alloys, and 
differences in the 800-hour alloy compositions in two of the alloys, 
did not facilitate accurate predictions beyond 500 hours. 

2. The Y-iayec widths and the weight of A1 consumed were likewise pre- 
dicted with good accuracy through 500 hours but with less accuracy 
beyond 500 hours. 

3. The experimental A1 concentration/distance profiles showed more 
curvature than those predicted for several of the alloys. The 
increased curvature could be the result of cross-term effects or. 
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more likely, a concentration dependence of the A1 diffusivity. 

4. Breakaway oxidation was correctly predicted to occur in the alloy 
with the greatest A1 content before an alloy with a lower A1 content. 



CONCLUSIONS 


Several conclusions may be drawn from the present study of the A1 trans- 
port during the cyclic oxidation of y+B NiCrAlZr(Y) alloys. The conclusions 
from the experimental results are: 

1. Oxide spallation results in an accelerated loss of A1 from the sample. 

2. Repeated oxide spallation results in a decreasing A1 concentration at 
the oxide/y-layer interface. 

3. The rate of A1 decrease at the oxide/y-layer interface is in propor- 
tion to the severity of the oxide spallation. 

4. Breakaway oxidation occurs when the A1 concentration at the oxide/ 
y-layer interface is approximately zero. 

5. The y-layer’ width thickens at an accelerated, non-parabolic rate. 

6. The loss of A1 from finite samples leads to diffusion path transla- 
tion, changes in the composition of both phases in the y+B region, 
and a related decrease in the volume fraction of the B phase. 

The conclusions from the numerical modeling study are; 

1. The diffusional transport of A1 in the y-phase layer can be modeled 
with acceptable accuracy to predict the y-layer thickness, the 
concentration/distance profiles, and the weight of A1 consumed. 

2, Breakaway oxidation can be predicted from the bulk alloy composition, 
approximate diffusion coefficients, and the requisite flux. 
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APPENDIX A 

TECHNIQUE OF QUANTITATIVE ELECTRON MICROSCOPY 

The Ni, Cr, and A1 concentrations in each of the NiCrAlZr(Y) alloys were 

determined by use of a scanning electron microscope (SEM), an energy-dispersive 

x-ray spectrometer (EDS), and a numerical correction scheme. The SEM was a 

JEOL JSM 35C operated at an accelerating potential of 20 KeV and a beam current 

-9 

of approximately 1 x 10 . amperes. The instrument was equipped with a beam 
current stabilizer. The x-ray source size, or approximate volume producing 
x-rays (65) was calculated to be 1.0- 1.5 ym. The EDS system consisted of an 
EDAX Si (Li) counting system coupled to a 711 series multichannel analyzer. 

The EDS system was connected to a Texas Instrument 733 ASR high-speed data 
terminal interfaced to a Data General N0VA2 mini -computer. The mini -computer 
contained both peak stripping routines and various correctional schemes. 

The correctional scheme used in this study accounts for the three main 
sources of interference between elements in an alloy. The three interference 
effects accounted for are: 

1. The variation in the generation of x-rays due to changes in the matrix 
composition and the incident electron energy. This effect is commonly 
referred to as the "atomic number effect" (Z) because of its depend- 
ence on the average atomic number of the matrix. 

2. The absorption (A) of the x-rays characteristic of one element by the 
other elements in the matrix. The degree of absorption is dependent 
on the depth distribution of the generated x-rays and on the absorp- 
tion coefficients of each of the other elements in the matrix. 

3. The secondary fluorescence (F) of x-rays generated by one element when 
it absorbs the primary x-rays generated by a second element. 

The equations to account for the above interference effects are commonly 
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referred to as ZAP corrections. The correctional scheme used in this study 
(66) was a modified version of the original computer program, designated FRAME, 
incorporating the ZAP corrections. Input to the program consisted of the x-ray 
intensity ratios (described later) for each element in the matrix, the char- 
acteristic x-ray line being used (K^ in the present study), the operating 
voltage, and the angles between the incident electron beam, sample surface, 
and the EDS collector. 

The typical sequence of steps to determine the concentration at a single 
point in the matrix is given below. The SEM was operated in the spot mode for 
a set time, usually between four and five minutes. The continuous spectrum 
with the characteristic peaks for Ni , Cr, and A1 was collected and stored by 
the EDS and multichannel analyzer. The continuous spectrum from the multi- 
channel analyzer was entered into the mini -computer, and the peaks of Ni, 

Cr, and A1 were stripped from the continuous spectrum by numerical routines 
stored in the mini-computer. The stripped peak intensities were divided by 
their respective pure element (Ni, Cr, and A1 ) intensities, measured under the 
same operating conditions. The resulting quotients, commonly referred to as 
the intensity, or k, ratios for each element were input back into the mini- 
computer along with the operating voltage of the SEM, the characteristic x-ray 
line used to determine the intensity rati os j and the appropriate angles as 
previously described. The computer program applied the ZAP corrections to the 
intensity ratios and listed the weight fraction for each element. In almost 
all cases, the sum of the weight fractions of Ni , Cr, and A1 was in the range 
0.95 - 1.00, most often between 0.98- 1.00. 

Nine NiCrAl alloys of unknown composition were analyzed by the above 
procedure to determine the accuracy of the ZAP correctional routines. Seven 
Y-phase alloys and two 3-phase alloys were fabricated and homogenized at 1200°C 
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for a week. The samples were sectioned and a portion of each sample was 
analyzed by atomic absorption at the NASA Lewis Research Center. The remaining 
portion of each alloy was analyzed following the procedure as described above. 
The compositions as determined by atomic absorption are compared with the 
compositions determined with the ZAP correction routines in Table A-1. As 
shown, the relative errors are the greatest for the A1 concentration; but the 
maximum relative error only amounts to less than 0.5 wt.% difference in the A1 
concentration. Therefore, the accuracy of the A1 concentration as determined 
in this study, at the worst, could be set at approximately ± 10% relative, 
which amounts to an absolute error of less than 0.5 wt.%Al 1.0 at.%Al). 


■ 4 » 

The position of the nine standards in relation to the measured phase 
boundaries of the NiCrAl system at 1200°C are shown in the phase diagram of 
Figure B-1 in Appendix B. 



TABLE A-1 


COMPARISON OF THE NI, Cr, AND A1 CONCENTRATIONS AS MEASURED 
BY ATOMIC ABSORPTION AND TECHNIQUE USED IN THIS STUDY * 

Atomc Absorption ZAF Correction % Relative AC 

Routines Error 


Ni 

Cr 


Ni 

Cr 

A1 

Ni 

Cr 


Ni 

9L 

M. 

Y phase 











62.57 

34.34 

3.09 

64.56 

34.59 

3.16 

3.2 

0,7 

2.3 

1.99 

0.25 

0.07 

66.54 

29.50 

3.96 

67.89 

29.20 

4,39 

2.0 

1.0 

10.8 

1.35 

0.30 

0.43 

70.28 

24.74 

4.98 

70.16 

24.67 

4.92 

0.2 

0.3 

1.2 

0.12 

0.07 

0.06 

76.41 

17.95 

5.64 

75,63 

17.62 

5.91 

1.0 

1.8 

4.8 

0.78 

0.33 

0.27 

73.40 

24.26 

2.34 

76.20 

24.50 

2.35 

3.8 

1.0 

0.4 

2.80 

0.24 

0.01 

79.83 

17.88 

2.29 

80.81 

17.54 

2.22 

1.2 

1.9 

3.0 

0.98 

0.34 

0.07 

90.34 7.10 

3 phase 

2.56 

89.30 

6.67 

2.33 

1.2 

6.0 

9.0 

1.04 

0.43 

0.23 

66.77 

14.84 

18.39 

65.92 

15.40 

18.58 

1.3 

3.8 

1.0 

0.85 

0.56 

0.19 

73.21 

8.00 

18.79 

73.49 

7.85 

18.53 

0.4 

1.9 

1.4 

0.28 

0.15 

0.26 


•k 

all concentrations in wt.% 

k k 

% relative error = (At. Abs, - ZAF)/ At. Abs.x 100 
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APPENDIX B 

EFFECTIVE NiCrAl PHASE DIAGRAM AT 1200°C 

An effective phase diagram for the NiCrAl system at 1200°C was determined 
and utilized throughout this study. Quantitative electron microscopy (Appendix 
A) in the y+3 region of various alloys yielded the concentration of they phase 
at the y/y + 6 phase boundary and the concentration of the B phase at the y/y+B 
phase boundary. The concentration of the y phase at the y/y+B phase boundary 
demonstrated a good fit to a single line when plotted on the NiCrAl phase dia- 
gram, whereas the concentration of the B phase at the B/y+B phase boundary 
showed more scatter about a single line. The experimentally determined data 
points and the effective phase boundaries, indicated as solid lines, are shov;n 
in Figure B-1. The concentrations and heat treatments are given in Table B-1. 

Two points of particular interest in this study can be noted by reference 
to the effective y/y + B phase boundary. First, the boundary shows no par- 
ticular dependence on the Zr or Y content of the alloys. In addition, the 
Ni-13.9Cr-25.lAl alloy, containing no Zr, shows the same approximate concen- 
tration at the phase boundary as the Ni-13Cr-25Al-0.05Zr alloy. Hence, there 
is no significant dependence of the y/y+B phase boundary on the Zr or Y con- 
tent of the examined alloys. Second, the concentrations of the y phase at the 
y/y+B phase boundary show the same fit to a single line regardless of whether 
the alloy was isothermally annealed as part of a diffusion couple, isothermally 
oxidized, or cyclically oxidized (where diffusion path translation caused a 
change in the concentration of the y and B phases with time). Therefore, it 
appears that the y and B phases were at thermodynamic equilibrium during the 
elevated temperature period of cyclic oxidation. 

The y/y+B and the B/y+B phase boundaries determined in this study are 
at lower A1 concentrations than the published NiCrAl phase diagram for 1150°C. 
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TABLE B-1 

ALLOY COMPOSITION AT THE y/Y+g AND B/y+B BOUNDARIES 


Alloy if 

Alloy 

Composition (at. 

%) 

Treatment 

1 

Cr 

13.3 

M 

24.7 

Zr (Y) 
0.05 

cyclically oxidized 100, 200 

2 

13.8 

25.0 

0.12Y 

500, & 800 hours, 
cyclically oxidized 100, 

3 

17.1 

28.4 

0.18 

200, 500, 1000 hours, 
cyclically oxidized 100, 

4 

19.0 

23.6 

0.04 

200, 344, 400, 500, 800 
hours . 

cyclically oxidized 100, 

5 

14.4 

23.7 

0.07 

200, 500, 800, 1600 hours, 
isothermally oxidized 200 

6 

14.4 

23.7 

0.07 

and 1000 hours . 
cyclically oxidized 800 

7 

13.9 

25.7 

0.0 

hours . 

isothermally oxidized 200 

8 

9.7 

17.2 

0.17 

hours . 

cyclically oxidized 200 

9 

14.4 

23.7 

0.07 

hours . 

isothermal anneal in excess 

10 

16.0 

17.5 

0.11 

of 100 hours (diffusion couple) 
isothermal anneal in excess of 

11 

20.8 

16.5 

0.02 

100 hours (diffusion couple) . 
isothermal anneal in excess of 

Standards 

Is 

35.9 

6.2 


100 hours (diffusion couple) . 
isothermal anneal for 1 week 

2s 

30.7 

7.9 

— , 


3s 

25.6 

9.9 

— 


4s 

18.6 

11.3 

— 



5s 


25.9 


4.8 
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TABLE B-1 (cont.) 


Alloy # Alloy Composition (at. %) 



Cr 

A1 

Zr(Y) 

6s 

19.2 

4.7 

— 

7s 

7.7 

5.4 

— 

8s 

13,6 

32.4 

— 

9 s 

7.3 

33.2 

— 
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The 1150°C phase diagram reported by Taylor and Floyd (67) is shov;n as dashed 
lines in Figure B-1. The differences between the measured y/y + B boundary and 
that reported are greater than the maximum ± 1.0 at.%Al accuracy of the con- 
centration-measuring technique (Appendix A). Taylor and Floyd also indicated 
that raising the temperature results in an increase in the A1 concentration of 
the y/y + 3 phase boundary. Hence, the difference in the y/y+B phase boundaries 
cannot be explained by reference to the 50°C temperature differential. A more 
accurate reassessment is necessary to clarify the discrepancy between the 
observed and reported phase boundaries. 
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APPENDIX C 

DIFFUSION COEFFICIENTS 


Introduction 

Diffusion phenomena in ternary alloys may be significantly different 
from those observed in binary alloys. In both alloys, atoms diffuse down 
their own chemical potential gradients, but whereas in the binary case the 
chemical potential gradient of a component increases with increasing concen- 
tration of that component, no such generalization can be stated for the 
ternary case. As a result, "uphill" diffusion, component diffusion up its 
own concentration gradient, or diffusion in a two-phase region is possible 
in ternary alloys. 

In a ternary system, four diffusion coefficients are required to describe 
the interdiffusion of atoms. The coefficients are of two types (64): direct 
diffusion coefficients which relate the flux of a component to its own con- 
centration gradient, and cross-term diffusion coefficients which relate the 
flux of a component to the concentration gradients of the other components in 
the system. It can easily be shown that two components are sufficient to 
analytically describe ternary diffusion (63). The appropriate diffusion 
equations relating concentrations, positions, and time for non-steady-state 
diffusion in a single phase is Fick's second law. The appropriate ternary 
form is: 


D 

dC, d ^11 
1 _ dx 


d 


dC^ 

dx 


dt 


dx 


dx 


dC2 

dt 


D 

d ‘^21 ^ 


dx 


D ^^2 
d ^22 


dx 


(C-la) 


dx 


+ 


dx 


(C-2a) 
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where: C^. = concentration of either component 1 or component 2 (i = 1,2), 

t = time, 

*^11’ ^22 ~ direct diffusion coefficients, 

^12’ ^21 ~ indirect diffusion coefficients, 

X = distance parameter. 

The terms containing 0^2 ^nd are commonly referred to as the cross terms. 
Pick's first law relates the flux of either component to the concentration 
gradients of both components. The appropriate form for Pick's first law is: 



dCj 

dx 

dCp 


11 

“^12 "dF 

(C-2a) 


dC, 

dC. 


21 

i 

dx 

^22 dx 

(C-2b) 


where 0^ and are the fluxes of component 1 and component 2, respectively. 
The other symbols are as described previously. Por the case of a moving 
boundary between a single-phase (y) and a two-phase (y+3) region, the 
appropriate mass balance at the interface is: 





for component 1 


(C-3a) 


(C^ - Cp ^ ^ for component 2 (C-3b) 

where: ^ = velocity of the y/y+3 interface, 

jT = flux of component i evaluated at the interface as approached 
from the y phase (i = 1,2), 

c9 = concentration of component i in the y+B region (i = 1,2), 

cT = concentration of component i at the interface as approached 
from the y phase (i = 1,2). 

A schematic of the concentration/distance profiles for the case where a y+3 
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NiCrAl alloy is diffusion bonded to a y-phase NiCr alloy is shown in Figure 
Cl. 

The effect of the cross terms on the diffusion behavior of NiCrAl alloys 
can be estimated by a judicious selection of y/y+B diffusion couples. It vas 
the intent of this particular investigation to estimate the effect of the 
cross terms on the diffusion of A1 and to investigate the possibility of using 
a single diffusion coefficient to describe the A1 transport. Towards this 
purpose extensive use was made of diffusion couple measurements on similar 
y+3 NiCrAl alloys by S. R. Levine of the Lewis Research Center (9). 

Experimental Procedure 

Various y/y+3 diffusion couples were utilized to investigate the dif- 
fusion behavior of three y + 3 NiCrAl alloys. The composition of the three 
y+3 alloys corresponded to the basic alloy compositions (minus Zr and Y) 
studied in the cyclic oxidation portion of this investigation. The y + 3 
alloys were cut, polished to 600 grit with SiC paper, and placed adjacent to 
various NiCr y-phase alloys in a molybdenum (Mo) canister (Figure C2). The 
canister and contents were annealed at 1200°C in an inert atmosphere for 
various times. The lower thermal expansion of the Mo canister caused a com- 
pressive stress to be placed on the couples resulting in diffusion bonding 
at the elevated temperature. The composition of the y+3 NiCrAl alloys and 
the y-phase NiCr alloys are listed in Table Cl. 

Diffusion of A1 from the Al-rich y+3 alloy into the NiCr alloy resulted 
in 3 dissolution and the consequent growth of a y-phase layer into the original 

The various y/y + 3 couples will hereafter be referred to by the short 
designation given in Table Cl. In this designation of i/j, i will refer to 
the y+3 alloy or "source" (high A1 ) , and j will refer to the y-phase alloy 
or "sink" (low A1 ) . 



o 


+ 


c. 

Figure C 


Distance 

■1 Y/Y+B Diffusion Couples 

a. NiCrAl source alloy and two possible NiCr sink 
compositions, b. schematic A1 concentration profile 
and, c. Cr concentration profile after annealing 
the diffusion couple. 



Figure C-2 Cross sectional diagram of Mo canister and diffusion couples 
used in determining the approximate diffusion coefficients. 
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TABLE C-1 


ALLOY COMPOSITIONS USED 

IN y/y+B diffusion couples 


Alloy Designation 

Ni 

Cr 

A1 

Zr 

Source 





1 

61.9 

14.4 

23.7 

— 

2 

57.9 

19.0 

23.6 

— 

3 

54.6 

17.1 

28.2 

— 

* 

4 

63.3 

12.7 

24.0 

— 

* 

4-Zr 

63.1 

12.7 

24.0 

0.16 

Sink 





W 

89.0 

11.0 

— 

— - 

X 

80.0 

20.0 

— 

— 

Y 

70.0 

30.0 

— 

. — 

Z 

60.0 

40.0 

— 

— 


* 

Reference 69. 
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two-phase alloy. A typical microstructure after a diffusion anneal is shown 
in Figure C-3. The diffusion couples were sectioned perpendicular to the 
original y/y+ 3 interface so that the Y-layer width could be measured. The 
couples were polished, etched, and the Y-layer width was measured by optical 
microscopy. The width was observed to increase at a parabolic rate in agree- 
ment with the previous study (9). The time dependence of the y-layer width 
for both studies is shown in Figure C-4. The A1 and Cr concentrations were 
measured in the y phase at the y/y+B interface in Figure C-1) for 

the 1/Y and 2/Y couples by the technique described in Appendix A. The limited 

data indicate that the measured concentrations (cX , ) are the concentra- 

Cr A1 

tionsatthe end of tie -lines passing through the y+B alloy compositions (see 
dotted line in Figure C-la). This behavior will be referred to as tie-line 
diffusion. Some Al^O^ precipitates appeared at the y/y+B interface, probably 
as a result of oxygen remaining in the sealed Mo canister. The appearance of 
the precipitate seemed to have little effect on the diffusion kinetics. 


Experimental Results and Discussion 

The approximately equal growth rates of the y-phase layer in couples 2/W. 
2/X, 2/Y, and 2/Z (2/_ couples) indicate the Cr concentration gradient had 
little effect on the A1 diffusion; i.e., is much less than Thi^ 

can be understood by reference to Eqs. C-2a and C-3a where component 1 repre- 


sents A1 and component 2 represents Cr. In Eq. C-3a, the y/y+B interface 
velocity (^) is approximately the same for the 2/_ couples. C^-j and C^-j 
can be assumed constant for the same source alloy (Figure C-la) . Therefore, 


from Eq. C-3a, the A1 flux ^ must be approximately equal for each of 
the 2/_ couples. The A1 concentration gradient (dC^-j/dx) can also be assumed 
to be essentially constant since the A1 concentration at the y/y+B interface 
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Figure C43 Diffusion couple 2/X after 132 hour anneal 
at 1200° C. The original YA + 3 interface 
position runs through the dark, colinear 
oxide precipitates. (150x) 
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Figure C-4 6 recession for the y+3 source alloys with various y 
sink alloys. Dashed lines indicate k values reported 
for couples 1/W, 2/W, and 3/W in Reference Q . Open 
symbols taken from Reference § . 
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Consideration of couples 2/X and 2/Z by similar reasoning implies the cross 

terms to be either dominant or significant with respect to the diffusive flux 

of Cr. The data in Table C-2 show that at a low A1 concentration, D„ and 

CrAl 

DcrCr similar in magnitude. 

Small additions of Zr have only a slight effect on the growth rate of 
the Y layer (69). Fi gure C-4 il 1 ustrates the slightly larger y-layer grov;th 
rate of couple 4/Y over the couple 4-Zr/Y (where the source alloy contains 
approximately 0.16 at.%Zr). Both couples have approximately the same y-layer 
growth rates as the couples with similar source alloys 1/Y, 1/W, and the 2/ 


couples. Significant in the 4/_ couples is the Kirkendall-type porosity in 
the y layer. Two points can be noted in regard to the porosity. First, 
there is roughly twice the amount of porosity in the y layer in the 4/Y 
couple where no Zr is present. Second, the location of the porosity is 
always on the y+B side of the original y/y+B interface. Micrographs illus- 
trating these two points are shown in Figure C-5. 

It may be inferred that the intrinsic diffusivity of A1 is greater than 
the intrinsic diffusivities of Ni or Cr in y+B NiCrAl alloys. Interdiffusion 
of Ni , Cr, and A1 occurs by the vacancy mechanism. In the four y + B alloys 
studied, A1 diffuses from the high-Al y+B source alloy into the MiCr sink 
alloys while the Ni and Cr diffusion is in the opposite direction (Figures 
C-lb and C-lc). The appearance of the porosity in the y + B side of the dif- 
fusion couple indicates a very rapid diffusion of A1 to the NiCr alloy re- 
sulting in a non -equilibrium concentration of vacancies on the y + B side. 

It is possible that the porosity is a result of vacancy coalescence into 
larger pores. This suggestion is supported by the fact that the intrinsic 
diffusivity of A1 is greater than the intrinsic diffusivity of Ni in Ni solid 
solution binary alloys (70,71). 
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Figtire C-5 Diffusion couples 4/Y ^a.) and 4-Zr/Y (b.) after a 
96 hour anneal at 1200 C. Couple 4/Y contains 
approximately twice the porosity of couple 4-Zr/Y, 
The original y/y+B interface is faintly observable 
in the left half of each micrograph, (200x) 
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The lower fraction of porosity in the Zr-containing alloy can also be 
explained by reference to a non-equilibrium concentration of vacancies in the 
Y + 3 side of the diffusion couple. Zr atoms are relatively large and would 
introduce compressive stresses about the Zr atoms in the Ni solid solution 
matrix. Formation of Zr- vacancy complexes could reduce the stress about each 
Zr atom. Hence, the formation of Zr-vacancy complexes would reduce the over- 
abundance of free vacancies, thus reducing the porosity in the Zr-containing 
couple. 


Discussion of Analysis 

It was decided to employ a binary analysis to determine a single diffusion 
coefficient approximating This decision was based on the previous 

observation that the A1 diffusion was most strongly influenced by the A1 con- 
centration gradient. An analytical solution was derived to describe the dif- 
fusion of A1 and the resulting parabolic growth of the y layer in y/y+ 3 
diffusion couples. Five assumptions were necessary for this analysis: 

1. The diffusion coefficients were concentration independent. 

2. The y/t + 3 interface was planar. 

3. No diffusion occurred in the y+ 3 region. 

4. The partial molar volume of A1 was independent of concentration and 
phase. 

5. The driving force for diffusion of A1 was solely due to the A1 concen- 
tration gradient; i.e., cross-term effects were considered negligible. 

The analytical solution (72) consists of an equation describing the A1 
concentration in the y layer and original NiCr alloy as: 


El + erf(;^)] 


^Ar^^ '"Al [1 + erf(Ti)] 


(C-4) 
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where: concentration measured from the original y/y + 3 interface, 

n = parameter specific to each y+ 3 alloy defined below, 

D = diffusion coefficient which approximates 
The Y-phase layer width is described by 


C = 2n/T5t (C-5) 

where ^ is the y-phase layer width. The n parameter is determined from a 
mass balance equation at the y/y+3 interface; 


n/i (C°^ - ^ exp(-n^) 

[1 + erf(n)] 

where all symbols are as previously defined. Eq. C-6 can be numerically 
solved for the n parameter for each of the y+3 alloys. From Eq. C-5 and the 
data in Figure C-4, the diffusion coefficient D may be determined. The 
parameter and diffusion coefficients are listed in Table C-3. 

The assumptions made to facilitate a binary analysis were considered 
acceptable for approximating the diffusion coefficient. The y/t + 3 interface 
was observed to thicken macroscopically as a planar interface, even though it 
was obviously not planar on the microscopic scale of the individual 3-phase 
particles. No concentration gradients were measurable in the y+3 region nor 
was there a detectable decrease in the volume fraction of the 3 phase near 
the y/y + 3 interface, both observations indicating insignificant diffusion 
in the y+ 3 region. The effect of a variable partial molar volume of A1 on 
the diffusion coefficient should be small and within the accuracy of the other 
measurements. The validity of ignoring the cross-term effects on the A1 dif- 
fusion has previously been discussed. No general statement can be made con- 
cerning the concentration dependence of the diffusion coefficient. It can be 
observed that for the three different alloys examined, the measured diffusion 



TABLE C-3 


INTERFACE 

COMPOSITIONS 

AND DIFFUSION 

COEFFICIENTS MEASURED ON THE Y/Y+3 

DIFFUSION COUPLES 

Couples 

(at. %) 

C^l (at. %) 

T) /j_ /T ^ \ 

7 C /t (10 cm /s) 

-10 2 

D (10 cm /s) 

1/W 


12.2-12.7^ 

0.38-0.39 

7.6 

12.5-13.2 

1/Y 

24.7 

12.2-12.7 

0.38-0.39 

7.2 

11.8-12.5 

2/W 


11.0-11.2^ 

0.32 

7.0 

17.1 

2/X 


11.0-11.2^ 

0.32 

7.6 

18.6 

2/Y 

29.6 

11.0-11.2 

0.32 

7.1 

17.3 

2/Z 


11.0-11.2^ 

0.32 

7.6 

18.6 

3/W 

35.0 

8.9-9. 4^ 

''0.2-0.22 

1.5 

7. 7-9. 4 


Assumed equal to 1/Y couple 

2 

Assumed equal to 2/Y couple 

Estimated from Ni-17Cr-28Al-0. 18Zr alloy examined in the oxidation section of this thesis 
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coefficients vary by less than a factor of three, not a very large difference 
with regard to diffusion coefficients, suggesting little concentration depend- 
ence over the range of concentrations examined. 

It is somewhat informative to compare the measured diffusion coefficients 
with the interdiffusion coefficients, D, in the Ni solid solution regions of 
the Ni-Al and Ni-Cr binary systems. The measured diffusion coefficients were 
in the range of the Ni-Al system shown in Table C-4. This observa- 

tion implies (64) that the chemical potential of A1 in NiCrAl alloys is 
principally a function of A1 concentration and virtually independent of Cr 
concentration. This suggestion is consistent with the previous observation 
that significantly less than Further support of the weak Cr- 

concentration influence is that the phase boundary between the y and the y + 3 
regions on the NiCrAl phase diagram lies almost parallel to the NiCr side, 
indicating a weak dependence of the chemical potential of A1 on the Cr concen- 
tration (64). 

The exact diffusion coefficient used in the computer model developed and 
used in this report is shown in Table C-5. The diffusion coefficients used 
for the Ni-17Cr-28Al-0. 18Zr and the Ni -14Cr-25Al -0. 12Y alloys were slightly 
less than the diffusion coefficients measured on the equivalent alloy not 
containing Zr or Y. Justification for this slight decrease was suggested 
previously in the section discussing the effect of Zr on the diffusion of A1 . 



p<^ 
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TABLE C-U 

COMPARISON OF Ni-Al AND Ni-Cr INTERDIFFUSION COEFFICIENTS AT 1200°C 


Ni-Al 


( 10 s ) 


Reference 


7.0 - 23,7 

71 

7.6 

68 

lU - 2h 

70 

,-Cr 

3.3 - ^.5 

37 

1.8 - 3.9 

73 

7.7 - 18.6 

This Study 


TABLE C-3 


DIFFUSION COEFFICIENTS USED IN THE NUMERICAL MODEL 


Alloy 


D (lQ~^^cm^/s) 


Ni-13Cr-25Al-0.05Zr 12.3 
Ni-lUCr-25Al-0 . 12Y 10.0 
Ni-lTCr-24Al-0.04Zr IT.l 
Ni-lTCr-28Al-0,l8Zr 6.2 
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APPENDIX D 
THE NUMERICAL MODEL 

A numerical model utilizing finite-difference techniques was developed 
to simulate the A1 diffusion during the cyclic oxidation of y+ 6 NiCrAl alloys. 
Finite-difference techniques have been well established as an effective tool 
in modeling complex diffusion behavior in both binary and ternary alloys. 

The numerical model developed in this study incorporates the complex boundary 
conditions at both the oxide/metal interface and the y/y+S interface observed 
to occur during the cyclic oxidation of y+ 3 NiCrAl alloys. 

Seven assumptions were made in the development of the numerical model. 

They were: 

1. The diffusion coefficients were concentration independent. 

2. The y/y + B interface was planar. 

3. The oxide/Y-layer interface was planar. 

4. The partial molar volume of A1 was independent of concentration and 
phase. 

5. The sample heat-up and cool-down times for each cycle were rapid; 
i.e., diffusion only occurred at 1200°C. 

6. Only Al^O^ was formed. 

7. The driving force for diffusion of A1 was solely due to the A1 
concentration gradient. 

The validity of each assumption is discussed in the numerical modeling section 
of this report. 

The selective oxidation of A1 in y''’B NiCrAl alloys results in the growth 


of a Y layer in the near-surface region. A solution of Pick's second law 
yields the A1 concentration/distance profiles in this layer. A unique solution 
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to Pick's second law requires knowledge of the boundary conditions at the 
oxide/y-layer interface and at the y/y+3 interface. A mass balance at each 
interface yields their position and velocity. The finite-difference equivalent 
of'Fick's second law, the mass balance at the oxide/y-layer interface, and the 
mass balance at the y/y+3 interface are shown in Table D-1. 

The finite-difference technique requires definition of a grid of points 
across the y layer. The grid points are equally spaced across the layer, each 
point being associated with a specific A1 concentration. Six to eight grid 
points were used in the present study. A linear concentration profile across 
a Y layer of 0.1- 1.0 pm width was assigned as an initial starting condition. 
The initial concentration profile and y-layer width were chosen so as to have 
negligible effect on the final solution. An initial concentration profile for 
six grid points is shown in Figure D-1. The grid spacing, Ax, is determined by 
dividing the y-layei" width - r,^) by n, the number of grid spacings desired. 
Interfacial A1 concentrations, and (at the oxide/y-layer interface and 
y/y + 6 interface) were assigned as the boundary conditions. 

The basic operation of the model is a sequential iteration of the finite- 
difference calculations for successive increments of time. The equations were 
therefore incorporated into a computer program to facilitate the repetitive 
calculations. Three fundamental calculations were performed in modeling the 
A1 diffusion during the cyclic oxidation of NiCrAl alloys. The first calcula- 
tion during the i'th iteration can be divided into two steps. First, the. A1 
concentration at the oxide/y-layer interface (C.^) and the A1 concentration at 
the y/y+ 3 interface were assigned for the particular oxidation time 

lAt. Second, the finite-difference form of Pick's second law operated on 
each concentration/grid point in the y layer, excluding the two interface 
concentrations previously assigned. The first fundamental calculation 
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TABLE D-1 

EQUATIONS USED IN NUMERICAL MODEL 


Partial Differential Equation 


Finite Difference Equivalent 


a. Pick’s Second Law 
D 


d^C 


dt dx 


(C*-C.) -2C.+C. , 

At Ax 


b. Mass Balance at T/Y+3 Interface 


(C -C ) 




0 ''yB' ■ At 


1 _ 


Ax 


c. Mass Balance at the oxide/Y layer Interface 




A? 


dt 


V J, 


2 - 


At 


= V D 


C -C 
^ 1 

Ax 


d. Murray-Landis Transformation 


none 


** 


O'? 


■■= c. + 
D 


S+1 ^:i-i 

2Ax 


j-1 Ag^ +n-j+1 A^2 
n n 


j=2,n 



refers to the concentration at the jth grid point at time t+At 
refers to the concentration at the jth grid point at time t after 


the M-L shift. 



21 


/ Phase 
Layer 


r + P 
Region 



Figure D-1 Initial Concentration Profile 
for 6 Grid Points 

n = 5 

• - C , at time = 0 
i 6 

- 0.1 - 1 . Oym 

= 0.0 

Ax = - ^ 2 ^/n = 0.02ym 

Cq = bulk composition of alloy 
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assigning the interface concentrations and operating on each grid point at 
time iAt is shown schematically in Figure D-2a, 

The second calculation during the same iteration, iAt, can also be 
divided into two parts. First, the flux of A1 from the y+3 region into the 
Y layer (J^ is calculated and the y/y+3 interface position is increased 
accordingly. Second, the flux of A1 passing from the y layer to the oxide 
(J 2 ) is calculated and the oxide/y-layer interface position (^ 2 ) '’s increased 
accordingly. The second fundamental calculation is schematically shown in 
Figure D-2b. 

The third calculation determines the new grid position based on the 
increased y-^aye^ width and adjusts the concentrations from the old to the 
new grid positions. This shift is sometimes referred to in the literature as 
the Murray-Landis (M-L) transformation (74). The third fundamental calcula- 
tion is schematically shown in Figure D-2c. The finite-difference form of 
the M-L transformation is listed in Table D-1. Following each iteration, the 
time is incremented to (i+l)At, the A1 concentration at both interfaces 
(C^ and at'e recalculated, and the three fundamental calculations are 

repeated. 

The total weight of A1 consumed in the oxidation reaction can be deter- 
mined from the flux of A1 entering the oxide at the oxide/y-layer interface. 
The A1 consumed during each iteration can be estimated as the product of the 
flux (J 2 ) and the time interval (At). The total weight of A1 consumed at a 
specific time is simply the sum of the weights of A1 consumed during each 
iteration up to the specified time. A conversion of the units of the weight 
of A1 consumed from that involving atomic percent to the more appropriate 
units of milligrams required the density of the y+ 3 alloys. A simple rule 
of mixtures of the density of each pure metal weighted by the atomic fraction 



Figure D-2a Main Iteration Scheme 

Operation of Pick's Second Law for the ith iteration 
A - Cg at time iAt 

A - C j concentration at grid j after calculation of 
Fick's Second Law 
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of each component in the alloy was used to estimate the density. Measured and 
estimated densities for some y and y + 6 alloys are listed in Table D-2. 

Recession of the metal surface at the oxide/metal interface is dependent 
on the partial molar volume of A1 in the y phase (see Table D-lc). The par- 
tial molar volume of A1 is a concentration-dependent variable. To determine 
the motion of the oxide/metal interface, the partial molar volume of A1 (V^-j) 
was estimated from the lattice parameter data for a Ni-20Cr-5Al alloy (67). 

3 

The determined value of 7.1 cm /mole can be compared with the molar volumes of 

3 3 3 

pure Ni , Cr, and A1 of 6.6 cm /mole, 7.2 cm /mole, and 10.0 cm /mole, respec- 

3 

tively. The partial molar volume was converted from units of cm /mole to 
units of inverse atomic percent to be compatible with the concentration units 
in the computer program. The calculation to determine the partial molar volume 
of A1 is shown in the following section. 

The partial molar volume of A1 determined for a Ni-20Cr-5Al alloy was 
taken as an average value for the y layer. Accounting for the motion of the 
oxide/metal interface results in an increase in the width of the y layer 
(measured from the original metal surface) and a corresponding increase in 
the weight of A1 consumed. The y layer increased only 6-12 percent in the 
four alloys studied when the oxide/metal interface motion was taken into 
account. For a specific alloy composition, increasing or decreasing the value 
of by 50% (more than spanning the range of the pure element molar volumes) 
resulted in less than a 5% change in the y-layer width. Hence, using an 
average value for the partial molar volume of A1 in the y layer was considered 
as an acceptable approximation in the present modeling study. 

The validity of the finite-difference technique was verified by comparison 


^(wt.%) 
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TABLE D-2 

, lEASURED AND ESTIMATED DENSITIES OF SELECT NiCrAl ALLOY 


Alloy 
(at. %) 

Ni-13C.r-25Al-0.05Zr 

Ni-19Cr-24Al-0 . 04Zr 

Ni-21Cr-16Al 

Ni-36Cr-6Al 

Ni-26Cr-10Al ' 

Ni-26Cr-6Al 

Ni-8Cr-8Al 


Measured Density. 
(gms/cc) 

7.4 

7.1 

7.4 
7.8 
7.8 

8.1 


Estimated Density 
(gms/cc) 

7.2 

7.1 
7.5 
7.9 
7.8 

8.2 


8.4 


8,4 
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with a closed-form analytical solution. Derivation of a closed-form solution 
required the same assumptions previously given for, the finite-difference solu- 
tion.' The general form of the analytical solution is given in Reference 72. 

In addition, the closed-form solution also required the assumption of a time- 
dependent concentration at the oxide/metal interface and at the y/y+B inter- 
face. As the number of grids used in the finite-difference solution was 
increased, the solution "converged" to the closed-form solution. The relative 
difference between the predicted interface position and the predicted weight 
of metal consumed, for' both solutions, decreased as the number of grid points 
used in the numerical solution was increased (Figure D-3). For the case of 
four or more grid points, the relative difference between the closed-form and 
finite-difference sol utions was less than 1%. 

The computation tim? for the numerical solution increased as the number 
of grid points in the y layer was increased. The time interval for each 
iteration (At) decreased as the number of grid points was increased, while 
the number of calculations during each iteration also increased. The number 
of necessary iterations for the indicated number of grid points is shown in 
parentheses in the lower plot in Figure D-3. Six to eight grid points were 
chosen as a balance between the necessary computation time and an acceptable 
error in the finite-difference model. 

Closely associated with a proper choice of the number of grid points is 
an appropriate choice of the time interval (At) during each iteration. The 
conservative condition 

At z: 0-25 Ax^ 

^ D 

was used for most of the modeling of the four alloys in this study. Doubling 
the time interval resulted in less than a 1% relative difference between 




of Grids 

. model to the analytical solution 
i number of iterations for the 
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solutions, whereas tripling the time interval produced no solution (similar 
to the case of one grid point in the upper plot in Figure D-3). 

A complete listing of the computer program incorporating the finite- 
difference solution is included at the end of this appendix. 

Partial Molar Volume of A1 

The partial molar volume of A1 at the concentration Ni-20Cr-5Al (vrt.%) 
was estimated by using a second-order Taylor series expansion. The partial 
molar volume of A1 at the concentration Ni-20Cr-5Al is determined by projecting 
a plane tangent to the three-dimensional volume surface at the concentration 
Ni-20Cr-5Al. The intersection of the tangent plant with the pure A1 axis is 
the partial molar volume of A1 at the concentration Ni-20Cr-5Al. The appro- 
priate form of the Taylor series expansion is: 



= V 



( 0.0 - 0 . 2 ) 



(0.0 - 0.75) 


lOOAl 'Ni-20Cr-5Al 


Ni-20Cr-5Al 


Ni-20Cr-5Al 


where: = intersection of the tangent plane with the pure A1 axis; 

this value is the partial molar volume of A1 at the concentra- 
tion Ni-20Cr-5Al. 

V = molar volume of a Ni-20Cr-5Al alloy. 

= change in the molar volume with changing Cr concentration 
evaluated at the concentration Ni-20Cr-5Al. 

= change in the molar volume with changing Ni concentration ■ 
evaluated at the concentration Ni-20Cr-5Al. 

The volume changes were estimated from the lattice parameter data of Reference 
67. The partial molar volume of A1 was initially determined in units of 
(A) /unit cell and was converted to cm /mole by reference to the fact of four 
atoms per unit cell in the face-centered cubic lattice structure and Avogadro's 
number of atoms in a mole of metal. 
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MAS Aa ALFLUX3 VSRSIOM 3 

NUMERICAL SOLUTION TO THE DlfFUSlONAL TRANSPORT OF AL IN 
NI"CR"AL ALLOYS UN0E8G02NS CYCLIC OXIOATION* 

THIS VERSION Of ALFLUX3 ALLOWS VARIABLE INPUT OF C2^ THE . 
AL CONCENTRATION AT THE OXIOE/GAMHA INTERFACE* 

PARABOLIC INTERFACE MOTION DOES NOT RESULT* 

PROBLEM; AL DIFFUSING OUT OF A TWO-PHASi GAMA-^-BETA 
REGION/ THROUGH A ONS“PHASE GAMMA REGION TO FORM 
AN OXIDE SCALE ON THE SURFACE, THE AL CONCENTRATION 
IN THE GAMMA PHASE VARIES FROM ALSURF AT THE OXIDE-GAMMA 
INTERFACE TO ALAB^ THE SOLUBILITY LIMIT OF AL IN THE GAMMA 
PHASE. A PLANAR INTERFACE IS ASSUMED WITH TWO MOVING 
BOUNDARIES# THE 6 AMMA/ GAMMA+BETA CXI1) AND THE 
OXIOE/GAMMA CXI2) INTERFACE 

THIS VERSION ALLOWS VARIABLE INPUT OF OXIDATION TIMES 


INPUT” 


42 

43 

44 

1 

OUTPUT “ 

C 

XI1 - INTERFAC 


C 

XI2 - INTERFAC 


C 

WM - HEIGHT OF 

47 

c 

FLUX OF AL AT 

48 

c 

AL CONCENTRATI 

49 

c 


50 

c 



C8#AL - INITIAL CR AND AL CONCENTRATION IN THE ALLOY CAT%) 
FORMAT » CR AT*% CCC 1“S5# AL AT«% CCC 7-11) 
XX,XX#XX.XX 

CRA8#ALAB - CR#AL SOLUBILITY IN GAMMA PHASE 
FORMAT “ SAME AS CR#AL ABOVE 
XX«KX#XX«XX 

ALSURF “ AL CONCENTRATION AT THE OKIDE/GAMMA INTERFACE 
FORMAT “ AL AT.5S CCC 1-5) 

XX. KX 

0 - DIFFUSION COEFFICIENT IN GAMMA PHASE 

FORMAT - CM*CM/SEC SCIENTIFIC NOTATION 
TP. 

TEMP “ OXIDATION TEMPERATURE IN DEGREES CENTIGRADE 
FORMAT - XXXX, IE, 1200 

V8A8 - PARTIAL MOLAR VOLUME OF AL IN THE GAMMA PHASE 
FORMAT “ XX, XX CCC/MOLE AL) 

NUM - THE NUMBER OF OXIDATION TIMES REQUESTED 
FORMAT - XX 

QXTIMECK) - THE SPECIFIC OXIDATION TIMES IN HOURS 
XXKXX.# XXXXX.# XXX.# XX.# XX, 

§ OF FIELDS EACH SEPcRATED BY COMMAS 


SET INITIAL CONDITIONS 


IMPLICIT REAL*8CA-H#L#M#0“Z) 

DOUBLE PRECISION J1 # J 2#NI #N I AB# J W# J GMS 

DIMENSION CC30)#CNEWC3C)#DI ST C30) # OXTIME C30) # CW C 30) 

READ *# CR#AL 
READ *# CRA3#ALAB 
READ *# CC15 
READ *# D 
READ *# TEMP 
READ *# VBAR 
READ *# NUM 

READ COXTIMECNN)# NN“1#NUfO 
READ *# N 
DO 15 NN=1#NUM 

15 0XTIM|<NN)=0XTIME<NN)*36D3 
NI=100-CR-AL 
N1AB=100-CRAB-ALA8 
XI1=10-6 

XI2=0.0 

DELT=,25*XI1*XI1F CD*N*N) 

VBAR = VBAR*7,5/CALAB*27,-«-CRAB*S2,+NIAB*58,7) 
0ENS=CNI*8,9+CR*7.2'FAL*2»7)/100 

MSSAL=0,0 

PRINT *#N 

is is it isifk itk is “kit k k k •k’k -kit is is k isis irk ^ is "kk k kirk 

PRINT HEADINGS AND INPUT SAMPLE 
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***************** iSrS:* ******** ********************* ******* ******* 

yRITE(6#150? CR^AL 

^ FORMATCJMI^ZX^J jtESULTS FOR SAMPLE NI-* / F5. 2/* ATX CR-%F5.2^ 

yRlTECdj'ISS) CRAB^ALAB 

F0RMAT(1H0<-5X,‘ CR/AL CONCENTRATION AT SAMMA/SAMM A+BETA % 

I , , • INTERFACE= *^2F6.2) 

WRITE(6/157) CC1J 

F0RHAT(6X/.‘ AL CONCENTRATION AT THE OXIDE/GAMMA INTERFACE* * /■ 

I F5» 2) 

WRITE(6^-160) TEMP 

F ORH AT <1 HO/ 5X/* OXIDATION TEMPERATURE COES C) » %FS,0) 

WRITEC6/200) D 

FORMATClHOfSX/* DIFFUSION COEFFICIENT CCM*CM/S£0 »%E14.4) 
yRITEC6/210) VBAR 

F ORMAT (1 HO/ 5X/* PARTIAL MOLAR VOLUME (1/ATX) *’/E14.4) 

******************************************************** 

SET INITIAL AL CONCENTRATION PROFILE 

USE LINEAR CURVE FOR INITIAL CONCENTRATION 

******************************************************** 

CCN+1 )=ALAB 
DELX=XI1/N 
A*(ALAB-C(1 ))/N 

DO 20 J=2/N 

ccj)«ca)-FCj-i)*A 

TIME=D£LT 


********************************************************* 
MAIN ITERATION LOOP FOLLOWS 

********************************************************* 

********************************************************* 
CALCULATE THE FLUX AND THE CHANGE IN XI/ DELXI 
********************************************************* 


50 J1 = CC<N“1)"4*CCN)+3*C(N+1 ))/(2*0ELX)*0 
0£LXI1=DELT*J1/(AL-C(N+1)) 
XI1=XI1+DELXI1 
IF<TIME.6T.3700CQ.) THEN 


Jtrkixncaui ajrijybUa/ i n cn 

CC1)*10.79-DSQRT<TIME/3600)*.012 

END IF 

J2=“ALPHA*D*<“C(3)+4.*C<2)-3.*C<1))/(2.*0ELX) 

DELXI2=»V9AR*J2*DELT 

XI2SXI2+0ELXI2 

0ELX=(XI1 "XI2)/N 

ALPHA=1«/(1.-VBAR*C(1)) 

AT.X TO 6MS/CC CONVERSION 

DO 55 11^1/3 

55 CW(II) = C(II)*27./ CC (II)-*27+NIAB*58.7+(100.-NIAB-CClI))*52) 
JW=D*(-CWC3)+4.*CWC2)“3«*CW(1))/C2.*DELX) 
ALAVG=CALAB+C(1))/2. 

CRAVG=100.”ALAVG-NIAB 

DENSG=CalAVG*2.7+CRAVG*7.2+NIAB*8.9)/100. 

JSMS=JW*DENS6 


ife * * * ^iir * * iSr ^ * A * i!r * * * *8: Sr iS- ^ ‘A A il^ * ^ * * * * * * iSr * * ifr * * * * * * ifr * 

SET OELT TO PROPER INCREMENT 

0ELT=»25*0ELX*0ELX/D 

TIME=TIME+DELT 

******** -Sr* ***it****il^*Sr*** *'*********** 

CALCULATE THE FLUX AND THE MASS OF ALUMINUM LEAVING IHE 
SYSTEM AT THE GAMMA-OXIDE INTERFACE 

************************************************ *********** 

MGFLUX=JGMS*1000. 

MGAL*MGAL'VABSCMGFLUX)*DELT 

******************************************************** 

USE THE MURRAY-LANOIS TRANSFORMATION TO SHIFT THE 
CONCENTRATIONS BACK TO THE CALCULATED CURVE 
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DO 60 J«2/N 

60 CNEW(J) = C(J)+CC<J+1 )-C (J-1) )/C2*DELX) 

1 *<DELXI1*<J-1 )/N+DELXl2*(N-J+1)/N) 

DO 65 J=2^N 
65 CCJ)sCNEW(J) 

********************************************************* 
CALCULATE NEW CONCENTRATIONS USING PICK'S 2ND LAW 
********************************************************* 

DO 70 J = 2<-N 

70 CNEW(j5=C(J)+<0*DELT*(C(J+1)-2*C(J)+C(J-1))/C0ELX*DELX)) 

DO 75 J=2,N 
75 CU) = CNEWCJ) 

IFCTIME.EQ. 360000) THEN 

CCN+1 >=0.0500* (TINE/ 3600) **0.5+9. 58 

END IF 

********* max iteration check ***************************** 

IFU.EQ. 200000) THEN , , , 

WRITEC6#100) 

100 FORMATC U^'MAX # OF ITERATIONS EXCEEDED') 


TIME=TIME/360D 
PRINT *^TIHE^0ELT^XI1 
STOP 

END IF 

1 = 1+1 

******** IF TINE EQUALS OXTIMECk) ************************* 

******** print results and stop *********** 

IFCTIHE.GT, COXTIME(K)“OXTIME(K)/200)) THEN 
DISTC1 )=XI2 
MM=HGAL 

TIMEH*TIME/3600 
, WRIT£(6/175) TIMEH 

175 FORMATC///5X, 'OXIDATION TIME (HOURS) =%F7.1) 

WRITE(6/’250) 

250 FORMAT(////5X/' CONCENTRATION PROFILE FOR GAMMA PHASE',- 
1 ///<,5X/'DISTANCE (C M) ' ,-1 OX,- ' AL (ATZ) 'r//) 

DO 90 J=2/N+1 
90 DIST(J)=<J-1)*DELX+XI2 
DO 95 J=1xN+1 

95 WRITE(6,3C0) DIST(J),-CU) 

300 F0RMATC8X, F9.5,r12X,-F6.3) 

WRITE(6,-400) WM , . . 

400 F0RMAT(lH0/5X,-'*= WM - AL LOSS FROM SAMPLE (MG/CM+CM) =', 

1 Ell. 5) 

WRITE(6/350) XII 

350 F0RMAT(1H0r5X/'XI1 POSITION (CM) ='rEl1.5) 

375 FORHATt^HQy5xJ^Xl2 POSITION (CM) ='^E11.5) 

WRITE(6,-450) MGFLUX 

450 F0RMAT(1H0,5X/' FLUX AT XI2 (OXIDE/GAMMA PHASE INTERFACE (MG/' 
1 ‘ CM*CM*SEC) ='/E9,3#///) 

PRINT *,l 

IF(K.GE.NUM) STOP 

K«K+1 

GO TO 50 

ELSE 

GO TO 50 
END IF 
END 
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